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Abstract of Dissertation Presented to the Graduate School 
of the University of Florida in Partial Fulfillment of the 
Requirements for the Degree of Doctor of Philosophy 

UF-IHEPA 00-01 

EXCLUSIVE MEASUREMENTS IN S ^ D*NNX 

By 

Antonio I. Rubier a 
August 2000 

Chairman: J. Yelton 
Major Department: Physics 

We report the first observation of exclusive decays of the type B — > D*NNX, 
where is a nucleon. Using a sample of 9.7 x 10^ BB pairs collected with the 
CLEO detector operating at the Cornell Electron Storage Ring, we measure the 
branching fractions B (5° ^ D*~ppn+) = (6.5ll:i ± 1.0) x 10-^ and B{B^^ 
D*^pn ) = (14. silo i 2.7) x 10~^. The charge conjugate process is imphed in 
the reconstruction of i?"— > D*~ p p 7r~^. However, in the reconstruction of 5"— > 
D*" p n, only the mode with the antineutron is used in our measurement because 
neutrons do not have a distinctive annihilation signature. 

Antineutrons are identified by their annihilation in the Csl electromagnetic 
calorimeter. Since we are unable to isolate a sample of antineutrons in data, we 
use antiproton annihilation showers in a A ^ pir'^ sample to define the antineu- 
tron selection criteria. We find a discrepancy for antiproton annihilation showers 
between the Monte Carlo and data, which we assume affects antineutrons as well. 
We increase the raw yield for > D*' p n by 21% to correct for this discrepancy. 
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The possible contributions from — > D*~ with p n and — > D*~ 

D*'^ with ^ -Dj"7 and Df pn are ehminated from the analysis by rejecting 
events with 1.91 GeV < M^^n < 2.04 GeV for a loss of 9% in the reconstruction 
efficiency. We fail to find evidence for the decay D'^ p n. 

We search for possible contributions to the resonant substructure of > D*~ 
p n and > D*~ p p tt"*" due to a heavy charmed baryon decaying strongly to p 
D*~ for > D*~ p p 7r+ and n D*~ for > D*' p n, as well as a resonance of 
the virtual W decaying to pp7r+. We also study the possible effect of feed-down A 
baryon contributions to the background for both modes, as well as the 5°— > D*~ 
p p TT'^ signal. No conclusive evidence is found for a measurable contribution from 
the aforementioned contributions to the resonant substructure. 

Antineutrons are used for the first time in the exclusive reconstruction of a 
B meson. By finding conclusive evidence for the existence of decay modes of the 
type B — > DNNX, we challenge the assumption that the B — > Baryons rate is 
dominated by decays to charmed baryons. 
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CHAPTER 1 
INTRODUCTION 



The study of the elementary particles at a wide range of interaction energies 
has occupied scientists since the discovery of the electron in 1896 [1]. Particle 
physics has evolved from a field involved in the discovery of new particles to one 
devoted to their systematic study. A logical structure, currently explained by the 
Standard Model of elementary particles [2, 3, 4], has been unveiled. 

The Standard Model, however, offers an incomplete picture of some experi- 
mental high energy physics results. The results we describe are amongst these. 
Despite some weaknesses, the Standard Model has withstood intense experimental 
scrutiny, and while some of the equations are currently not calculable, evidence 
has not been found for physics beyond the Standard Model. Particles acquire their 
masses in the Standard Model via the Higgs mechanism, that requires the existence 
of a massive gauge boson, the Higgs boson, that has yet to be found. 

We first introduce particle physics phenomenology -particles and their decays, 
which is the information most often used by the practitioner of experimental high 
energy physics. A discussion of the current status of B physics phenomenology 
follows. Finally, we review B Baryons previous to our results, concluding with 
an overview of our work and its significance. 
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1.1 Matter 

Matter is composed of three types of indivisible constituents: leptons, gauge 
bosons, and quarks. Leptons and gauge bosons are found individually in Na- 
ture. Quarks combine in two types of arrangements to form hadrons. The first 
arrangement is of the form quark- antiquark, and is called a meson. The second 
arrangement, three quarks, is called a baryon. Mesons and baryons, collectively 
known as hadrons, comprise all the known forms of matter consisting of quarks. 

Hadronic matter is said to be colorless. A color is assigned to each quark in a 
meson or a baryon. The three quarks in a baryon each have a different color— red, 
green, blue— and the combination of all three colors yields a colorless hadron. The 
quark and the antiquark in a meson form a color- anticolor pair (e.g. quueQuue)- 

1.1.1 Hadrons 

There are three families of quarks, each consisting of two types of quark: an 
up-type quark, with electromagnetic charges +2/3 that of the electron, or (+2/3) 
Qe, and a down-type quark, with (-1/3) Qg. Every type of quark is called a flavor 
of quark. 

The first and lightest family consists of the up (u) and down (d) quarks. The 
next family, with heavier quarks, is the charm (c) and strange (s) family. Even 
heavier still is the third family: the top (t) and bottom (b) quarks. All are shown 
in Table 1.1. 

For every matter constituent there is an anti-matter constituent with opposite 
electromagnetic charge and equal mass, as shown in Table 1.2. 
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Table 1.1: Quark families 



up-type (+2/3) 


up: u 


charm: c 


top: t 


down- type (-1/3) 


down: d 


strange: s 


bottom: h 



Table 1.2: Anti-quark families 



up-type ( 


-2/3) ge 


u 


c 


t 


down-type 


(+1/3) qe 


d 


s 


b 



Quarks are not found individually in Nature. Their masses can be estimated 
by the mass spectrum of mesons and baryons that has been measured to date. 
In Table 1.3 we show current best estimates of the lower and upper limits for 
quark masses. The estimated masses for quarks are useful in current Standard 
Model calculations. However, the quark mass estimates have large uncertainties, 
especially in the case of the up and down quarks. The values quoted in this Table 
and the next values quoted in this section are from the 1998 Review of particle 
physics [5]. The mass of the top quark rrit = 173.8 ± 5.2 GeV. 



Table 1.3: Estimated quark masses 



Quark 


Low mass 


High mass 


u 


1.5 McV 


5 MeV 


d 


3 McV 


9 McV 


s 


60 MeV 


170 MeV 


c 


1.1 GeV 


1.4 GeV 


b 


4.1 GeV 


4.4 GeV 



The proton, for example, has the quark content (uud), and one of the pions, 
the TT'^, has the quark content {ud). Each type of quark is considered a distinct 
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flavor. The group of mesons containing a charm quark and a hght antiquark (one 
of s, 11, d) is called the charm mesons. 

Meson and baryon masses are known to varying degrees of accuracy, as shown 
in Table 1.4 for selected mesons, and Table 1.5 for selected baryons. 



Table 1.4: Selected meson masses 



Meson 


Quark content 


Mass 




ud 


139.56995 ± 0.00035 MeV 




su 


493.677 ±0.016 MeV 




cu 


1864.6 ± 0.5 MeV 


5° 


bd 


5279.2 ± 1.8 MeV 



Table 1.5: Selected baryon masses 



Meson 


Quark content 


Mass 


proton 


uud 


938.27231 ± 0.00028 MeV 


A 


sud 


1115.683 ±0.006 MeV 


Ac 


cud, 


2284.9 ± 0.6 MeV 




bud 


5624 ± 9 MeV 



1.1.2 Leptons. 

The leptons, the e, the fi, and the r, are fundamental particles. Each can be 
produced with an accompanying neutrino, vi, with I — e, /i, or r. Neutrino physics 
in the near future will yield a better understanding than is currently available. 
Neutrinos have been thought to be massless and not to mix (that is, each neutrino 
flavor was thought only to interact with its lepton partner) . The recent observation 
of neutrino mixing, however, implies that neutrinos have mass [6]. Unlike quarks. 
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leptons are observed as single particles with well-measured masses. In Table 1.6 
we show the current mass measurements. Large differences in masses are found 
between the electron, the muon, and the tau. 



Table 1.6: Lepton masses 



Lepton 


Mass 


electron: 
muon: /i^ 
tau: 


0.51099907 ± 0.00000015 MeV 
105.658389 ± 0.000034 MeV 
1777.05lg;l GeV 



Helicity is the orientation of a particle's momentum vector in respect to its 
spin. Helicity can be right-handed or left-handed for all fermions except neutrinos, 
which are only left-handed. Antineutrinos, likewise, are only right-handed. 

1.1.3 Gauge Bosons 

The mediators of the physical forces, the particles that allow decays to take 
place, are called gauge bosons. In Table 1.7 we outline their masses and the 
types of interaction that each mediates. I — e,ii,T, u— up-type quaxk{u, c,t), and 
d=down-type quark((i, s, h). 



Table 1.7: Some characteristics of bosons 



Boson 


Force 


Mass 


Decays 


photon 


electromagnetic 


massless 


7 ^ /+/- 


gluons 


strong 


massless 






weak 


80.41 ± 0.10 GeV 






weak 


91.187 ±0.007 GeV 


W+ l+l-,uP,q+q- 



1.1.4 Spin and Statistics 

Matter is also characterized by the statistics obeyed. Leptons and quarks are 
fermions, obeying Fermi statistics, an example of which is the Pauli exclusion 
principle for electrons occupying the same shell in an atom. The spin of leptons 
and quarks is (±|). Bosons obey Bose statistics, which allow an infinite number of 
particles to occupy the same energy state. Bosons have integral spin (0, ±1). The 
different relative alignments of the spins of the individual quarks, together with 
the addition of angular momentum, results in a large number of possible states. 

1.1.5 The CKM Matrix 

The interactions of quarks from different flavor families are suppressed with 
respect to those within the same family. In order to correct for this discrepancy, 
the characteristics of a given decay as prescribed by weak theory are adjusted using 
the flavor mixing 3x3 unitary and complex matrix V, the Cabibbo-Kobayashi- 
Maskawa matrix [5] , which transforms the weak d, s, b quark states to the physically 
measured d' ,s' , b' quark states, to arrive at a correct theoretical understanding of 
a weak decay for all quark flavors: 



Kd (0.9745 - 0.9760) K.(0.217 - 0.224) K6(0.004 - 0.013) ^ 



CKM 



V;d(0.217 - 0.224) V;,(0.9737 - 0.9753) ^6(0.036 - 0.042) 
\4d(0.004 - 0.013) \4,(0.035 - 0.042) 146(0.9991 -0.9994) 



with 
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(A 



s 



V 



s 





The diagonal elements of V are 1, implying that decays which involve these 
CKM matrix elements are Cabibbo-allowed, whereas all other decays, which involve 
off diagonal elements, are Cabibbo-suppressed. For example, a 6 — > cW~, W~ — > 



1.1.6 Symmetries 

A particle and its antiparticle are said to be symmetric under CPT transfor- 
mations, where C is charge conjugation, P is parity, and T is time reversal. The 
inversion of parity acts like a mirror in the inversion of space coordinates. CPT 
symmetry is valid for all forces. 

A subset of this symmetry is CP. The product of a charge conjugation and a 
parity inversion affect a particle the same as its antiparticle. When CP is violated, 
a particle will prefer a different subset of decays than its antiparticle. Such is the 
case for neutral kaons.The and states in weak theory are different from 
the physically observed strong states, the Ks, or X-short and the Kl, or X-long, 
each with its respective antiparticle. Short and long refer to the respective decay 
lifetimes. Charge symmetry is obeyed implicitly in neutral kaon decay, but parity 
is violated. CP violation is evident in the difference in decay rates of Ks^s and 
Kls to final states of two and three pions, the former being even under parity. 



cs decay has a much larger decay rate than a 6 ^ cW ,W cd decay. 
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and the latter odd under parity [7] . CP violation in B decays and will be the focus 
of many studies in the near future. 

1.2 Decays 

All hadrons, except the proton and electron, are unstable and decay to lighter 
hadrons, leptons, and bosons. The top quark decays well before the time required 
to form hadrons (baryons or mesons). Therefore, the bottom group of baryons 
(e.g. A;,, with quark content bud) and mesons (e.g. Be, with quark content be) is 
the group of hadrons with the largest masses currently found in Nature. 

1.2.1 Weak Decays 

The b quark in the B meson decays via the weak interaction. The b quark 
decay is often accompanied by the strong decays of soft gluons which allow for the 
hadronization of a large number of final states. The Fermi Electroweak theory, 
which aimed to explain neutron f3 decay, n pe~i7e, as shown by the Feynman 
diagram Figure 1.1, introduced the neutrino to particle physics, serving as the 
precursor of the Standard Model. A full theoretical understanding of this decay 
was accomplished by the Standard Model, in which this decay is mediated by 
the W vector boson, not present in the Fermi Electroweak theory. The decay 
n — > pe~z7g is correctly described as the quark level process d — > W~u, followed by 
W~ — > e^harvei as shown by the Feynman diagram in Figure 1.2. 




P 



Figure 1.1: A Feynman diagram for neutron beta decay in Fermi Weak Theory 



id 



neutron 




proton 



Figure 1.2: A Feynman diagram for neutron beta decay in the Standard Model 
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Examples of weak semi-leptonic decays of mesons, analogous to neutron (5 
decay, are K'^ n^e'^Ug, — > K'e^Ue, and S° D~e^Ve. The weak semi- 
leptonic decays of baryons follow a similar pattern. These decays are referred to as 
semi-leptonic decays since the final products are a combination of a weak leptonic 
decay and quark hadronization. Leptonic weak decays, in which the unstable 
particle annihilates into an pair, examples of which are tt"*" and 
K'^ jJi^Vfj, are the most theoretically well understood type of decays, as it lacks 
any final state hadronization. 

1.2.2 Strong Decays 

The weak decays of the heavy (charm and bottom) quarks in mesons and 
baryons often involve the secondary emission of strongly-interacting soft gluons. 
While we understand the weak component of these decays using the Standard 
Model, the strong component is not calculable. 

The strong interaction hadronization process is not well understood theoreti- 
cally when soft, or low momentum, gluons mediate the decay. The Standard Model 
is based on perturbation expansions which rely on convergence. A process involv- 
ing soft gluons yields equations that are no longer perturbatively convergent. This 
stumbling block has prevented us from understanding many details of unstable 
particle decay, particularly for the case of heavy hadrons, in which the large mass 
available in the decay implies a very large number of possible final decay products 
from an equally large number of soft gluons. 

The D* meson, for example, which we reconstruct in this analysis, is a spin 
1 meson that decays via the strong interaction. Two possible Feynman diagrams 
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for this decay are shown in Figures 1.3 and 1.4. Our ignorance about soft gluon 
strong interactions forbids us from knowing the proportion of the total decay rate 
due to any one diagram. 




d d 



Figure 1.3: A Feynman diagram for D* D^tTso/j 




Figure 1.4: A second Feynman diagram for D* D^ngoft 
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1.3 B Meson Decays 

The B meson, which is the topic of our study, can decay to a large number of 
hghter particles by various mechanisms, some of which are detailed in Table 1.8. 
The b — > cW~ type decays account for more than 95% of the B decays that are 
possible in the Standard Model. The combined semileptonic decay rate for the 
three lepton families is ^ 25%, with hadronic decays accounting for almost all of 
the remaining decay rate. 



Table 1.8: Some B meson decay mechanisms 



Quark-level mechanism 


Sample decay mode 


Semileptonic decay: 
b cW-, W- l^ui 

Hadronic decays: 
b^cW-,W- cs 
b^ cW-,W- ^ ud 


B Dl-ui 

B DDs 
B^ Dn 



The W boson is colorless. Its decay to quarks constrains the color of both 
quarks to cancel. The number of Feynman diagrams for a given B decay by the 
color of the quarks. Whereas decays mediated by an external W boson allow for 
any of the three possible colors (for example, the 7r~ in B~ — > D^Tr~, as shown in 
the Feynman diagram in Figure 1.5). Those in which the W boson decays internally 
limit the color of all the quarks to be the same as that of the parent meson, as 
shown in the Feynman diagram in Figure 1.6. The latter quality is referred to as 
color suppression. Decays that only have an internal W boson-mediated diagram, 
such as B^ D^tt", are color-suppressed decays. 
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The decays b — > (c, u)W'' are mediated by the W~ vector boson, while decays 
of type b — > (s, d)^ and b (s, d)g are mediated by neutral bosons. All decays 
except b — > cW~ contribute a small fraction of the total decay rate. However, it 
is for many of these rare decays that our current phenomenological understanding 
allows for the closest agreement between theory and experiment. 

It is theoretically allowed, and has been experimentally measured, that a 
meson can oscillate to a B° before decay, allowing for neutral T(45') — > BB events 
with either 2 S°'s, or 2 S°'s [8]. 

1.3.1 Quantum Chromo dynamics 

The Standard Model unifies the electromagnetic and the weak interactions. It 
also encompasses the strong interactions in the form of Quantum Chromodynam- 
ics (QCD) [9, 10]. The strong coupling constant is smallest for short range 
interactions, or for large momentum transfer, a quality of QCD called asymptotic 
freedom. 

The QCD Lagrangian is given by [11]: 



Qred 



^QCD = {qred, Qblue, Qgreen) {^1^ D ^ 



m) 



Qblue 



Qgreen 



where the covariant derivative is: 



'A 
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are the SU(3) flavor matrices, and are the eight color gauge fields. G^i, 
is the gluon field-strength tensor. 

1.3.2 Heavy Quark Effective Theory 

The properties of the D and B mesons, which are composed of a heavy quark 
and a light anti-quark, have been described using heavy-quark symmetry [12, 13] 
by Heavy Quark Effective Theory (HQET). In a meson, a quark and an antiquark 
are confined to a bound state in a cloud of virtual quarks and gluons which need 
to be incorporated into any calculation. In the case of a heavy meson, such as a S 
meson, the heavy quark has a substantially larger mass than the light antiquark. In 
HQET, the heavy quark is independent of the light anti-quark. HQET assumptions 
simplify the Standard Model equations, allowing, for instance, the comparison with 
theory of experimental values of Vub and Vcb- 

The effective Lagrangian that is used to characterize B decays is given by 
[14, 15]: 



where Gp-i the Fermi constant, is 1.17 GeV ^, and Jqq is the charged weak 
current given by: 



JCC = {^e, l^n, 1^t)Y 



( \ 



{uL,CL,tLh''VcKM 



Sl 



18 



The assumption that the heavy quark mass uiq is effectively infinite is used to 
simphfy the QCD Lagrangian. The heavy quark and the hght quark are decoupled, 
and the effect of the cloud of virtual light quarks, light antiquarks, and gluons is 
assumed to be small enough to be ignored. The QCD Lagrangian, of which C^ff 
is a simplified version, is simplified to: 

where D'^ is the QCD covariant derivative. 

In the limit uiq oo the strong interactions of the heavy quark become 
independent of its mass and spin and the effective Lagrangian is further simplified 
to: 

^HQET = hviv ■ Dhy 

where is the effective heavy quark field and v is the hadron's velocity, which 
is close to that of the heavy quark. 

In the calculation of HQET quantities, the strong interaction effects that are 
non-perturbative are grouped into a form factor that includes a dimensionless 
probability function, the Isgur-Wise function ^{v ■ v') [16], where v and v are 
respectively the initial and final velocities in a scattering or decay process. An 
example of the role this function plays is the elastic scattering of a S meson. The 
hadronic matrix element for this process is: 



:^{B{v)\h^,YK\B{v)) = i{v ■ v'){v + v'y 
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where and b^i are the heavy quark fields. The heavy quark symmetry used 
in B physics phenomenology represents substantial progress in the theoretical de- 
scription of B decays. In the next section we discuss semileptonic B decays, for 
which HQET has been successfully used to derive decay rates. 

1.3.3 Semileptonic Decays to Mesons 

The combination of large branching fractions and high reconstruction efficien- 
cies have allowed experiments such as CLEO to measure several semileptonic B 
decays with high accuracy [17]. This wealth of experimental results has allowed 
phenomenologists to compare theory and experiment. The decay kinematics of a 
specific B semileptonic decay dictate the type of form factor contributions to the 
decay rate. In the case of S — > D*W, for example, there are no ^ corrections to 
the decay rate. This behavior, which is explained by Luke's theorem [18], implies 
that the HQET-derived decay rate for B — > D*^ has low sensitivity to higher 
order perturbative corrections as well as non-preturbative effects. 

1.3.4 Hadronic Decays to Mesons 

Whereas HQET has been useful in describing semileptonic B decays, an equally 
accurate description of hadronic B decays using HQET has only recently begun to 
be pursued for two-body decays to mesons [19, 20, 21]. Whereas in semileptonic B 
decays one of the two currents in a matrix element is weak, and therefore calculable 
to all orders in perturbation theory, in the hadronic case we have matrix elements 
of four-quark operators with hadronic uncertainties due to the exchange of gluons 
and quarks. 
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In energetic hadronic two-body B decays hadronization is assumed to take place 
after the quarks that form each of the two hadrons have traveled sufficiently long 
distances for there to have been no significant exchange of gluons or quarks between 
them. This decay characterisitic is referred to as factorization, in which the matrix 
elements of four-quark operators factorized to independent current elements for 
each hadron. By using the operator product expansion (OPE) [22, 23], the weak 
interaction effects are treated as separate from the long range strong interaction 
effects. The HQET effective weak hamiltonian for h ^ c,u transitions after this 
procedure is given by: 

^e// = ^(K6[ci(//)gi'^ + C2(At)g2i + h.C.) + Q;,^^ + Qpenguin 

where ci, 2(/i) are scale dependent Wilson coefficients. Qi,2*-i Qb^ui and Qpenguin 
are, respectively, four-quark operators for 6 ^ c, 6 ^ and penguin decays [20]. 
The relative strength of each type of operator as well as the Wilson coefficients are 
decay-dependent . 

Applying the factorization hypothesis to H-eff to, for example, the decay am- 
plitude of the decay — > D'^n', results in: 

Afactorized = ^KbK/ai (tT" | ( 0) (D+ 1 (c6),,,,,, | E^) 

where oi can be verified with experiment. As the number of hadronic two-body 
B decays and the accuracy with which their decay rates are measured increases 
in the near future, it will be possible to test the decay rates derived using HQET 
similarly to how it has been done for semipletonic B decays. 
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1.4 S — > Bciryons 

A distinctive feature of the B meson system is that the large mass of the b-quark 
allows for many of the weak decays of the B meson to include the creation of a 
baryon-antibaryon pair. The as yet unmeasured decay p n bars this feature 

from being unique. In Table 1.9 we outline the B — > Baryons decays allowed in 
the Standard Model. 



Table 1.9: Representative B — > Baryons decay mechanisms 



Quark-level mechanism 


Sample decay mode 


Semileptonic decays: 
b^cW-,W- l-ui 
h uW-, W- I'ui 
Hadronic decays: 
b cW-,W- cs 
b^ cW-,W- ^ cs 
b^ cW-,W- ^ ud 
b^ cW-,W- ^ ud 

b sg, 
b^uW-,W- ^ud 


B Kpl'vi 
B ppl'vi 

B — >• J/ipK, J/ifj ^ Baryons 

B SeAc 
B AcpTT 

B^^ D*- p n 
B — > Kp,B ppK 

B — >• pflTl 



1.4.1 Results to Date 

The B ^ pX was measured by CLEO to be 8.0 ± 0.5 ± 0.3%[24], assuming 
B — > pX — B ^ nX. Based on this number we expect roughly 8% of all B 
mesons in our data to be a S ^ Baryons event. The Ac, Sc and Sc charmed 
baryons have been inclusively measured in B decays [25, 26]. Upper limit exclusive 
measurements have been reported to date for B — > Acpl'vi with I — e, /i [27], and 
selected two-body rare decays [28] . 
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Based on the inclusive measurements reported to date, B — > Baryons decays 
to date have been expected to be produced predominantly in decays of the type 
B —>■ AcpX, and this is the type of decay that has been exclusively reconstructed 
to date [29]. A typical Feynman diagram for a 5 — > KcpX decay is shown in Figure 




Figure 1.7: A Feynman diagram for B — > Ac^pn . 



1.4.2 The Argument for B [D]NNX modes 

Combining the recently measured value B{A^ — pK~7i~^) = 5.0 ± 0.5 ± 1.2% 
[30] with estimates of the product branching fraction B{B AcX) x B{Ac — > 
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pK'n'^) — 0.18% [31] we can determine that B A^pX modes account for only 
3.6% oi B ^ Baryons, which is approximately half of the total B — > Baryons 
rate, as measured in the inclusive B — > pX measurement [24]. 

Based on our current knowledge of S — > Baryons, there must be processes 
other than Ac production that contribute to the B Baryons rate. It is our goal 
to find evidence for B — > Baryons decays which do not involve a Ac. Dunietz 
[32] suggested that modes of the type B — > DNNX, in which D represents any 
charmed meson, and N a proton or neutron, are likely to be sizeable. B — > DNNX 
final states can arise either from the hadronization of the W boson into a baryon- 
antibaryon pair, or from the production of a highly excited charmed baryon that 
decays strongly into a baryon plus a charmed meson. CLEO previously reported 
an inclusive upper limit for B{B — > DNNX) at the 90% confidence level of < 4.8 
%[24]. 

1.4.3 Thesis Overview 

In this thesis we will attempt the exclusive reconstruction of two specific B — > 
[D]NNX decay modes, >■ D*~ p p and B^^ D*^ p n. Typical Feynman 
diagrams for B^^ D*~ p n and D*~ p p tt"*" are shown in Figures 1.8 and 

1.9, respectively. 

The choice of these two modes is guided by the following criteria: 

1. D* mesons have the lowest signal contamination among the D mesons. 

2. Both decays can occur via external W decay. Although this characteristic is 
not a principle, to date only the b —> c decay modes that have been measured 
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share this characteristic. The reasons for the predominance of these decays 
are not known. 

3. These are the two modes with the lowest decay daughter multiphcity, which 
translates into the highest reconstruction efficiency. 

4. > D*~ p n and > D* p p have low combinatoric backgrounds. 

We report here, for the first time, evidence for decays of the type B 
[D]NNX, and present measurements of the branching fractions B{B^^ D*~ p 
p 7r+) and B{B^^ D*~ p n). The charge conjugate process is imphed in the recon- 
struction of B^^ D*^ p p TT"*". However, in the reconstruction of B^^ D*~ p n, 
only the mode with the antineutron is used in our measurement because neutrons 
do not have the distinctive annihilation signature. These measurements invalidate 
the previous assumption that B — > Baryons is dominated by Ac decays, while es- 
tablishing evidence for the existence of a new type of B decay mechanism with a 
sizeable decay rate. 

The thesis is divided as follows: 

In Chapter 2 we describe the CLEO detector. We place special emphasis on 
the electormagnetic calorimeter, which we use to select antineutron candidates. 

In Chapter 3 we outline the selection criteria for the B°^ D*~ p p tt"*" and 
D*~ p n decay daughters. We place special emphasis on our selection of 
antineutron showers due to the novelty of their use. 

In Chapter 4 we describe our measurement of B^^ D*~ p p 7r+ , in which we 
use a recontruction technique used to reconstruct decays in which the energy of all 
the decay daughters is well determined. 
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In Chapter 5 we describe our measurement of p n , in which we 

use a recontruction technique which is similar to that used to reconstruct other B 
decays with missing energy. In our case the missing energy is due to the antineutron 
candidate. 

In Chapter 6 we conclude by summarizing our results, stressing on their sig- 
nificance, and outlining possible decay modes that we believe are important and 
measurable with the expectedly larger datasets available to future studies. 





Figure 1.8: Two Feynman diagrams for D* p n 




Figure 1.9: A Feynman diagram for D* p p 



CHAPTER 2 
CLEO II DETECTOR 



The CLEO II detector [33] is centered at the interaction point resulting from 
the collision of an electron beam and a positron beam, each at ^ 5.3 GeV beam 
energy, located at the Cornell Electron Storage Ring (CESR). The beams collide 
almost head-on, resulting in a center-of-mass total energy of ^ 10.6 GeV. It is in 
the energy range near this value that the T resonances, with the quark content bb, 
are produced. As shown in Figure 2.1, the T(15') to T(45') resonances are found 
in the energy range 9.44 GeV to 10.62 GeV. The T(45') is above the threshold for 
strong decay to BB pairs to take place. Close to 100% of the T(45') decay rate is 
to nearly equal numbers of charged and neutral BB pairs. 

Data is taken at the T(4S') resonance (to be referred to as ON resonance) 
to study B decays. Because approximately 2/3 of the ON resonance data are 
composed of events in which the initial qq pair is not an T(4S'), a sample of data 
is also taken 60 MeV below the resonance (to be referred to as OFF resonance) in 
order to subtract the non-T(45'), or continuum component of the ON resonance 
data. The e~^e~ annihilation at or near the T(4S') resonance yields a wide variety 
of possible final states, some which are shown in Table 2.1. 

The decay product of an e'^e" annihilation is called an event. Because final 
states differ greatly in cross section (the frequency with which events of a given 
topology are produced) some event types are produced with high frequency, while 
other event types are produced with low frequency. Electronic triggering on an 
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Figure 2.1: Cross section into hadrons from the collision of e+e beams at CESR 
as measured by the CLEO II detector in the energy range 9.44 GeV to 10.62 GeV 
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Table 2.1: Some e^e annihilation final states 

e+e~ 77 
e+e~ (with l=e,/x, r) 

e^e" T{AS) — > BB (requires an energy threshold) 
e+e~ qq (followed by hadronization) 



event-by-event basis is used to select some or all of a given type of event. Triggering 
allows us to record only the events that we are interested in studying, most of which 
are of types produced with low cross sections. 

2.1 Sub-detector Components 

In order to get enough information out of an event, we use an ensemble of 
sub-detectors, each yielding incomplete information about the event. When the 
information from all sub-detectors is combined, we have sufficient information to 
measure useful physics properties. A front view and a side view of the CLEO II 
detector are shown in Figure 2.2 and Figure 2.3. 

From the innermost to the outermost (with respect to the beam pipe, which is 
located at the center of the detector) , the sub-detectors are: 

1. Vertex detector 

(a) PTL (precision tracking layers) detector, used during the earher part of 
data recording. These data will be referred to as CLEO II data. 

(b) 3- layer SVX (silicon vertex detector), used during the later part of data 
recording. These data will be referred to as CLEO II. 5 data. 




Figure 2.2: Front view of the CLEO II detector 
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Figure 2.3: Side view of the CLEO II detector 



33 



2. Drift chamber. 

3. TOF (Time-of-Flight) detector. 

4. Electromagnetic calorimeter. 

5. Muon detector. 

The volume including all but the muon detector is enclosed in a 1.5 Tesla 
superconducting magnet. An important feature of this magnet is the uniformity of 
its magnetic field, which ensures that charged particles bend uniformly regardless 
of where in the detector the particle travels. A clear introduction to detectors as 
well as experimental methods in high energy physics is found in Perkins [34]. 



2.2 Tracking System 

A charged particle traversing a magnetic field in the presence of charged wires 
in a chamber containing gas will ionize this gas as it loses energy. We measure 
the time at which this process takes place as well as the energy collected by each 
wire. These measurements allow us to know the position of the particle in time 
and the energy released at a number of points, or hits, along its trajectory. For a 
given momentum, the rate at which a particle loses energy along this trajectory, 
measured as dE/dx, is dependent on its mass, thus allowing us to separate protons, 
kaons, and pions. The TOF (Time-of-Flight) of a particle in a scintillating medium 
is also dependent on its mass and momentum. TOF measurements yield a second 
way to separate protons, kaons, and pions. 
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2.2.1 PTL Detector 

The PTL detector is an inner drift chamber composed of six layers of straw 
tubes. There are 64 axial wires for each layer, and there is a half cell stagger 
between sequential layers. The PTL detector does not measure the longitudinal, 
or z-axis, position of the particle. The PTL transverse position measurements are 
more precise than those from the drift chamber. 

2.2.2 SVX Detector 

In later running of CLEO II, the PTL drift chamber was replaced by a 3-layer 
SVX detector capable of longitudinal as well as axial measurements [35], each 
measurement taking place on the two sides of each of the sihcon wafers. The 
radii of the SVX layers are 2.35 cm, 3.25 cm, and 4.75 cm, for layers 1, 2, and 3, 
respectively. It is composed of 96 wafers arranged into 8 octants of 12 wafers each, 
with 26,208 data readout channels. The intrinsic resolution from e'^e~ — > 
events at normal incidence is 29 //m. 

The improved measurement resolution of the SVX detector in comparison with 
the PTL detector allows for more accurate determination of the event vertex. This 
advantage is most useful to lifetime studies, yet it does not affect greatly the results 
presented here. 

2.2.3 Drift Chamber 

The drift chamber system (the main drift chamber and the vertex detector), 
together with the SVX or PTL, are used to measure the momentum of charged 
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particles. Some vertex detector and drift chamber parameters are shown in Table 
2.2. The beam pipe is located at radius 3.5 cm in the CLEO II data and at radius 
2.0 cm in the CLEO II. 5 data. 



Table 2.2: Vertex detector and drift chamber parameters 



Detector 


Layers 


Radius (cm) 


Wires per layer 


PTL (CLEO II only) 


6 


4.7 to 7.2 


64 


Vertex Detector (VD) 


10 


8.4 to 16 


64 (first 5), 96 (second 5) 


Outer Drift chamber 


51 


17.5 to 95 


96 to 384 



The r — and z measurement resolutions for each of the sections are shown on 
Table 2.3. 

Table 2.3: Drift chamber resolution 



Detector 


r — (f) resolution 


z resolution 


PTL (CLEO II only) 
Vertex Detector (VD) 
outer drift chamber 


90 fj,m 
150 iim 
110 //m 


N/A 
0.75 mm 
3 cm 



The Vertex Detector (VD) is bounded by concentric inner and outer cathode 
strips which provide z measurements. The segmentation of the VD cathode strips 
is 5.85(6.85) mm along z, which is the beam direction, on the inner(outer) cath- 
ode. Segmented cathodes also bound layers 1 and 51 of the outer drift chamber. 
Segmentation is about 1 cm along z. 

2.2.4 Momentum and Anguleir Resolution 

There are two factors that affect the track momentum resolution: 
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1. The error in the measurement of the track curvature due to the hit-level 
measurement error in drift distance. This resolution component is parametrized 
by a term linear in pt, the transverse momentum. 

2. Multiple scattering at material boundaries which cause the track trajectory 
to deviate from a helix. This resolution component is parametrized by a constant 
term. 

The parametrization for CLEO II data is, in GeV: 

(|i)2 = (o.0011j9i)2 + (0.0067)2 

In Table 2.4 we show this resolution in MeV for selected values of pt- 

Table 2.4: Momentum resolution {5pt) for the CLEO II data at selected values of 
Pt 



Pl 


<^Pl 


0.5 GeV 
1.0 GeV 
1.5 GeV 
2.0 GeV 


3.3 MeV 
6.8 MeV 
10.4 MeV 
14.1 MeV 



The angular resolution is measured using e'^e~ — » events, in which the 

typical Pt = 5.0 GeV. 54> = 1 mrad and 59 = 4 mrad. We expect S4> and 69 to 
be higher for the tracks we use in our analysis since muons have a much lower 
probability of multiple scattering than do other charged particles. 

2.2.5 dE/dx Measurements 

dE/dx is a function of particle mass and momentum, since p = m/37, where 
P — v/c. The degree of separation we are able to achieve is shown in Figure 2.4 
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for the CLEO II data. The 51 layers of the outer, or main, drift chamber— to be 
referred to as DR— are used to measure the specific ionization energy loss (dE/dx) 
of particles. 

A mixture of argon-ethane gas was used in the main drift chamber for the 
CLEO II data. This mixture was changed to helium-propane during data taking 
for CLEO II. 5, which allowed for an improvement in the dE/dx resolution, resulting 
in better charged particle separation. In Figure 2.4 each of the particle bands has 
been plotted after a subtraction of some of the higher dE/dx data points in a given 
particle band. 

The main drift chamber has contiguous cells each with a sense wire surrounded 
by field wires, as shown in Figure 2.5. Overall, there are three field wires for every 
sense wire in the main drift chamber. A number of corrections are applied to 
optimize the resolution: 

1. Dip angle saturation: tracks perpendicular to the sense wires have the highest 
density of ionization along the z direction. The amount of collected charge 
is reduced by electric shielding for these tracks. 

2. Drift distance: varies depending on the field configuration of each cell. 

3. (r, 0) entrance angle: its magnitude as well as its sign is field dependent. 

4. Axial-stereo layer: cells for axial and stereo layers have a different field de- 
pendence. 
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Plot not available 

For information on the CLEO detector, please see: Y. Kubota et al., Nucl. Instrum. 
Meth.A320:66-113,1992 



Figure 2.4: dE/dx vs. track momentum 
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VD Outer Shell 



Figure 2.5: Main drift chamber (DR) wire arrangement 
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2.2.6 Time-of-Flight Measurements 

Time-of-flight detectors surround the DR. A bar of organic scintillating material 
5 cm thick has photomultiplier tubes at each end for the barrel detectors, and at 
one end for those in the endcaps. A time measurement is made, with a 154 ps 
resolution for the barrel, and 272 ps resolution for the endcaps. From the time and 
distance travelled, a quantity is defined. 1/(5 varies by particle type, as shown 
in Figure 2.6. 

2.3 Electromagnetic Calorimeter 

The electromagnetic calorimeter is used to measure the electromagnetic energy 
deposition of charged and neutral particles. It is composed of 7800 Csl crystals. A 
clustering algorithm is used to combine the energy deposition in a crystal region, 
which is called a shower. 

The information from detector components is used in our analysis in a way 
that is consistent with previous CLEO II analyses with the exception of measure- 
ments from the electromagnetic calorimeter. The electromagnetic calorimeter has 
been used previously to measure electron and photon energy deposition. In our 
measurement of > D*~ p n we are required to select showers that are consis- 
tent with being due to antineutrons annihilating with the Csl. The antineutron 
selection procedure is successful for the first time at CLEO. 
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Plot not available 

For information on the CLEO detector, please see: Y. Kubota et al., Nucl. Instrum. 
Meth.A320:66-113,1992 



Figure 2.6: Time-of-Flight vs. track momentum 
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2.3.1 Dimensions 

The calorimeter is witliin tlie 1.5 Tesia magnetic field. All crystal faces are 
at 1 m from the interaction point, facing it in such a way that showers reach all 
crystals at normal incidence. A partial diagram showing some of the barrel and 
one of the endcaps is shown in Figure 2.7. Each calorimeter crystal is 5-cm x 
5-cm X 30-cm, where the latter is the length of the crystal. The choice of thalium 
doped Csl for the calorimeter crystals took into consideration factors such as cost, 
resistance to cracking, high density, and short radiation length (1.83 cm). Because 
the calorimeter is 16.4 radiation lengths deep, ~ 1% of the energy of a 5 GeV 
electron leaks out of it. There are 6,144 barrel crystals and 828 crystals for each 
endcap. 




IP PTL 



Figure 2.7: Layout of CLEO II detector showing barrel and endcap calorimeter 
sections 



43 

The crystals are arranged in a rectilinear grid with care taken to have a min- 
imum of material between crystals. The material in front of the barrel section is 
described by radiation lengths in Table 2.5. 

Table 2.5: Radiation lengths of material before barrel section of the calorimeter 



Material 


Radiation length (%) 


beam line to layer 1 of DR 


2.5 


argon-ethane gas and rest of DR 


0.4 


DR outer cathode layer 


1.0 


outer DR wall 


2.0 


TOF counters 


12.0 



There is more material in front of the endcap calorimeter sections, which de- 
grades measurement quality. We use endcap information in the 7r° reconstruction, 
but not for the antineutron selection criteria. The central barrel, which covers 71% 
of the polar angle (45° < 9 < 135°), has less material in front of it, thus providing 
the measurements with highest resolution. 

2.3.2 Clustering 

The algorithm involved in associating a group of nearby cells with energy above 
a threshold in the calorimeter is called clustering. An average of 430 crystals have 
energy recorded in a hadronic event (either BB or continuum). The raw ADC 
count measurements need to be calibrated using a crystal-to-crystal calibration. 
The sample used is Bhabba events (e+e~ — > e+e~). This sample has high statistics 
and a beam energy constraint. The electrons and positrons in Bhabba events 
deposit almost 100% of their energy in the calorimeter. Crystal noise, which is 
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in the range of a few MeV, does not affect this sample appreciably. The set of 
constants obtained from the crystal-to-crystal calibration needs to be updated 
only every few months. The main change is the preponderance of a few crystals to 
become noisy. 

A clustering algorithm for CLEO should accomplish the following: 

1. Match tracks to showers. 

2. Allow 7r° reconstruction. 

3. Allow separation of photons and K^s from antineutrons. 

There are two clustering packages in the CLEO II data: CCFC, and XBAL. 
Both accomplish these three requirements. Even though we use CCFC, using 
XBAL would have yielded results similar at the few % level. Our discussion of 
clustering is hmited to CCFC. A discussion of clustering using XBAL would be 
similar to ours. 

In order to separate shower clusters from noise, only cells with energy > 10 
MeV are considered as candidates for the center of a shower. Only cells that are 
at most two cells away from a cell in the cluster are added to the cluster, which 
allows cells without an energy measurement to be located inside a cluster. The 
number of cells N that is considered a cluster is a function of the total energy of 
the cluster. For example, 4 cells correspond to a 25 MeV shower, 17 cells to a 4 
GeV shower. 

The position vector, which contains the directional cosines, is determined using 
a weighted sum of the energy measured with each cell, using the geometric center 
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of each to define a centroid. Two corrections, one lateral, and another longitudi- 
nal, are applied to this position vector. These corrections account for calorimeter 
segmentation as well as depth into the Csl crystal of the shower center. 
The cluster energy E and angular resolutions, using Monte Carlo, are: 

Barreh ^[%] = |^ + 1.9 - l.OE 
Endcap: ^[%] = -p| + 2.5 
Barrel: (7^[mra(i] = ^ + 1.9 
Endcap: (7^[mra(i] = + 7-3 
Barrel: (Jo[mrad] = O.SG^sin{9) 
Endcap: aQ[mrad\ — ^ + 5.6 

Using the example of a photon in the barrel, we hst in Table 2.6 energy and 
azimuthal angular resolutions at two cluster energies. Their low values are a testa- 
ment of our ability to properly reconstruct electromagnetic showers. A by-product 
of this ability is our success in separating baryon-antibaryon annihilation showers 
from electromagnetic showers. 

Table 2.6: Energy and angle resolutions for a photon in the barrel at two values 
of cluster energy 



Cluster energy 


F, 


[%] 




(T(i,[mra(J\ 


100 MeV 


3 


8 


11 


5 GeV 


1.5 


3 
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2.4 Muon Detector 

By the time a particle has traversed all the sub-detector components enclosed 
within the muon detector, all particles except muons have deposited most of their 
energy in the calorimeter and/or have been deflected by the magnetic field within 
the drift chamber. Since the muon detector is, above everything else, a large 
volume of iron, only muons are expected to pass through a significant section of 
the muon detector. 

The nuclear absorption length for a muon in iron is 16.8 cm. It is considerably 
less for all other charged tracks we study. The maximum depth of iron to be 
traversed by a particle is between 7.2Ai and lOAj. We do not use the muon detector 
in our analysis. 



CHAPTER 3 
PARTICLE SELECTION 



In order to reconstruct > D*" p p and > D*~ p n we need to 
have low background samples of all the decay daughters. Protons, antiprotons and 
pions are selected as single tracks, the mD*-m£)0 distribution is used to select D*'s, 
and shower selection criteria is used to select antineutrons. We select the B decay 
daughters from ON resonance hadronic events. The ON resonance candidates 
can be fully or partially reconstructed B mesons, or background candidates from 
ON resonance continuum events. The OFF resonance candidates can only be 
continuum events as the total event energy is below the energy threshold for BB 
production . We use the Mn_fit histogram package [36], and its function fitting 
utility, MINUIT [37], both of which are widely used by experimental high energy 
physicists. 

To suppress the already small number of continuum background candidates in 
our reconstructions, we select events using the parameter R2GL [38]. R2GL is 
a measure of how the momentum is distributed for the event. A high value of 
R2GL corresponds to continuum events, in which the initial quarks hadronize to 
form a two-body decay. The BB momentum distributions are the result of two 
separate two-body or higher decays that can be most approximately described by 
a spherically symmetric momentum distribution, which tends to yield a low value 
of R2GL. 



47 



48 

3.1 Data Sample 

The CLEO data is collected at the T(45') resonance (ON resonance) and 60 
MeV below the resonance (OFF resonance). Roughly 2/3 of the ON resonance 
data is continuum and the remaining 1/3 is BB pairs. We show in Table 3.1 the 
integrated luminosity by dataset. 



Table 3.1: The CLEO data 



Dataset 


ON luminosity 


OFF luminosity 


BB's 


CLEO 11 
CLEO II.5 


3.14 /r' 

6.03 fh-^ 


1.61 fh-' 
2.94 fb-^ 


3.30 X 10^ 
6.35 X 10^ 



3.2 Monte Carlo Sample 

We use Monte Carlo generated events to model the data in order to study the 
efficiency of our selection criteria as well as the effect of backgrounds. Two types 
of Monte Carlo samples are used: 

1. Generic Monte Carlo, in which many decays are generated using our current 
knowledge of B decays. 

2. Signal Monte Carlo, in which one of the B mesons in the event decays via a 
predetermined mode, such as > D*~ p n . The remaining B meson in the 
event decays according to the prescription used to generate generic Monte 
Carlo events. 

The generic Monte Carlo sample consists of B decays of type b — > cW~. A 
detailed explanation of the criteria used for the generic Monte Carlo sample is 
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found in the Appendix of [24]. This sample does not include any B — > [D]NNX 
decays because until now we have had no evidence for their existence. The B — > 
Baryons decays in the generic Monte Carlo sample are of the type B — > niy+, in 
which n — > Ac~N, where N is an antiproton or an antineutron. The conjugate 
modes are also generated. 

3.3 Track Selection 

A combination of measurements is used to select quality tracks that are most 
likely to be the particle type we need for a given reconstruction-pion, kaon, or 
proton. The tracks we use in our reconstruction are particles that are stable (in 
the case of protons) , or decay outside our detector (in the case of kaons and pions) . 
Therefore, a track is assumed to begin near the interaction point, which is within 
a solid volume defined by the uncertainty with which we can measure the beam 
spot. This volume is a few microns wide in x and y and a few hundred microns 
wide in z. 

3.3.1 Fitting Algorithm 

Tracks are fitted using DR hits by a Kalman algorithm that is based on the 
assumption that the track is a helix inside a vacuum. The change of trajectory 
due to material interaction is taken into consideration in the track fit. The fit is 
performed from both ends of the track. The inward fit begins a liit-to-hit swim, 
adding hits to form a track, analogous to the way beads are strung to form a 
necklace. Where there is more than one hit to a layer of the drift chamber, the 
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hit with the better resolution is chosen. The outward fit begins at the outermost 
layer of the drift chamber, performing a hit swim into the drift chamber. 

3.3.2 Drift Chamber Track Variables 

Variables we use which are defined with main drift chamber measurements 
(DR) are: 

1. PQCD: signed 3- momentum of track. This vector does not have a correction 
assuming the mass of a particle type. 

2. ZOCD: distance of closest approach along the z-axis. The z measurement res- 
olution is considerably worse than the r — measurement resolution. How- 
ever, ZOCD is effective in separating tracks from the interaction point from 
material interaction tracks and long decay particle daughter tracks. 

3. DBCD: radial distance of point of closest approach to the interaction point. 
DBCD is usually more accurately determined than ZOCD. 

All tracks are required to pass the following track quality cuts: 

1. PQCD < 1 GeV: DBCD (in meters) < 0.005-(0.0038) x \PQCD\. 

2. Low momentum tracks are more likely to have a poorly measured DBCD. 
PQCD > 1 GeV: DBCD (in meters) < 0.001. 

3. ZOCD < 0.05 meter. 
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3.3.3 The TRKMNG Package 

A software package is used to reject duplicated tracks from curlers [39] . A curler 
is a low momentum track which is bent more than one revolution inside the drift 
chamber. A curler spirals inside the drift chamber as the momentum of the particle 
decreases. The tracking algorithm will often assign a track number in the list of 
tracks for each half-revolution. The TRKMNG package selects the half-revolution 
that is most hkely the closest to where the track began to curl, rejecting all other 
half-revolutions. The variable used in the TRKMNG Package is tng(track) > 0. 
All negative numbers of tng(track) are for the redundant half-revolutions. 

3.4 Particle Separation 

By population size alone, all tracks are pions. The purpose of particle sepa- 
ration, or particle ID, is to separate tracks that are likely to be something other 
than a pion from the overwhelmingly more numerous pions. Protons are produced 
with a significantly reduced rate, therefore constituting a small background to 
all other charged particles. The predominant backgrounds to protons are kaons, 
pions, electrons and muons. For the case of kaons, pions represent the largest back- 
ground. The kaon in D^^K^t:~ has more phase space available than the kaon in 
D^^K^t:~t:^ and >ir+7r~7r"'"7r~, which allows its kinematic separation from 
backgrounds. Even for the low Xp in £)*'s from B decays, the kaon in D^^K^t:' 
can be easily separated from the pion by the decay kinematics. 

Charged particle identification is accomplished by combining the specific ion- 
ization (dE/dx) measurements from the drift chamber with time-of-fiight (TOP) 
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measurements. A normalized probability ratio Lj is used for charged particle sep- 
aration. 

Li is definded as Lj = Pi/{Ppion + Pka<m + Pproton), and Pi is the particle hy- 
pothesis Xi probabihty combining dE/dx and TOP measurements, with Xi defined 
as: 

- ,2 _ r {dE/ dx)measured- {dE / dx)expected ^ 2 , [ {TOP) 

measured {'^^■^^ expected 2 
^ C^dE/dx '■ C^TOF ^ 

o'dE/dx,TOF are the deviations from the mean and vary by particle type. We 
define the respective Xi for each particle type using the crdE/dx,TOF for the identical 
particle type. Lkaon, for example, is defined using the (TdE/dx,TOF^ for kaons. In 
Table 3.2 we show how our criteria varies, depending on how much information 
is available. By defaulting to a pion we are setting Pi = 1.0. Useful dE/dx 
information is available for the vast majority of tracks passing the track quality 
requirements outlined in section 3.3.2. Low momentum particles (< 200 MeV) curl 
before reaching the TOP detectors. 



Table 3.2: Particle identification criteria 



dE/dx 


(JtOF 


Action 


< 3 


< 3 


use both 


< 3 


> 3 


use adE/dx only 


> 3 


< 3 


use (Ttof only 


> 3 


> 3 


default to piou 



We place loose requirements of Lpion > 0.001 on all pion candidates, and the 

kaon candidates in D^^K^n^ . Kaon candidates in D'^^K^ti^h^ , and D^^K^Ti^Tr^ir 
are required to satisfy Lkaon > 0.4. Proton candidates in D*~ p p tt~^ and 
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> D*~ p n are required to satisfy Lproton > 0.9. In addition, electron back- 
grounds to protons and antiprotons are rejected using the variable cut R2ELEC 
< 0. R2ELEC is a logarithmic probability that a particle is an electron. Negative 
values of R2ELEC are least likely to be for an electron. R2ELEC is defined using ^ 
from the electromagnetic calorimeter, which is very close to 1 for electrons, which 
leave all of their energy in the calorimeter, as well as a non-negligible number of 
protons. 

3.5 7r° Reconstruction 

The TT^'s used to reconstruct D^^K^n'ii'^ decay via 7r° 77. The vr^'s in this 
decay are soft, with momenta < 500 MeV. The decay angle for 7r° 77 is large 
enough at this momentum range for the individual 7's not to overlap. On average 
there are 7 tt^'s in a BB event. The average number of tt'^'s, as of yet unmeasured, 
should be near this number. Since shower multiplicity increases with decreasing 
shower energy, as 7r° momentum decreases, the combinatoric shower background 
increases. A minimum momentum of 100 MeV is required of our tt" candidates to 
suppress this background. 

A shower is considered a 7 if the ratio between the summed energy of a 3 by 
3 grid of shower cells over the same are extended to a 5 by 5 grid of shower cells 
is almost 1. This ratio is referred to as E90E25. Most photons have most of their 
energy concentrated within a 3 by 3 grid of shower cells. 

We require at least one of the 7's to be in the good barrel section of the 
calorimeter, which is defined as \cos{0^)\ < 0.71. When one of the 7's is not in 
the barrel, we require that it not be from the overlap region between the barrel 
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and endcap sections of the calorimeter or near the beam pipe, both of which yield 
degraded shower energy measurements. We also require the candidate 7r° mass to 
be within 2.5 a (tt*^ ct ~ 5 MeV) from 135 MeV central value. 

3.6 D* Reconstruction 

Finding D* candidates with high efficiency and low backgrounds is an important 
part of this analysis. In the decay D* — >■ D^Tisoft mD*-TnDO is only slightly larger 
than the mass of the tt^. Since the mD'-fnoo resolution is smaller than the mo* 
resolution, we use the former to select D* candidates. We verify our D* yield with 
a previous analysis. 

In order to reconstruct B D*X with D* D^ngoft correctly, we check our 
yield for the scaled momentum Xp range 0.25 < Xp < 0.35. It is in this range that 
we find ~ 50 % of D*'s to be generated in signal B^-^ D*^ p p and > D*~ 
p n Monte Carlo. We check our yield with a previous CLEO result [40]. 

3.6.1 The KNLIB Fitting Package 

The KNLIB fitting package is a collection of kinematic routines used for track 
vertexing [41]. We use it to vertex D^'s and D*^s. D°'s are vertexed using the 
KNLIB routine KNVTX3. A new invariant mass is calculated with this vertex, and 
only the mass range 1.70 to 2.03 GeV is allowed to pass on to D* reconstruction. 
To form a D* vertex, we use KNLIB routine KNBVXK which combines the 
decay products and the tt^ track to make a vertex-constrained mass from which 
mD*-mDO is calculated. We use KNBVXK with the flag IVOPT set to 1 for the tt^- 
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i.e. the tt^ is constrained to pass through the vertex formed by the other particles, 
but not used to determine that vertex. The latter method offers a noticeable 
improvement over other options tested, including use of tt^ track when vertexing, 
as well as using the I?° vertex instead of the individual D° tracks when vertexing 
the D*. We place a loose cut of < 100 on the D* vertex fit for all three modes. 

3.6.2 Fit Optimization 

We reconstruct D°'s, selecting candidates near the mean mass. After the ad- 
dition of a TTs, we select candiates within a range of m£)*-m£)o values. We use the 
1998 Review of Particle Properties values for m£)0, 1.8646 GeV, and mD*-moo, 
0.1454 GeV, as our central values [5]. The double Gaussian fits model mD*-TnDO 
better than the single Gaussian fits. 

The optimization procedure is as follows: 

1. The background is fitted with a second order Chebyshev polynomial. 

2. The signal is fitted with a double Gaussian distribution. 

3. The continuum subtracted data moo distribution is fitted using abs(mD*- 
mDO-0.1454) < 0.002 GeV, which is a wide cut. 

4. The widths, as well as the means of each of the double Gaussians are allowed 
to float. The best fits are obtained with non-zero Am ban's. Am ban is the 
difference between the mean value of each of the two Gaussians. We integrate 
the double Gaussian, keeping 95% of the signal symmetrically from the PDG 
mean. 
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5. The final set of double Gaussian m^jo and mD*-mDO cuts from data are used 
in the B reconstructions, and the same cuts are applied to the Monte Carlo 
to find cmc-i the signal Monte Carlo reconstruction efficiency. 

Shown in Figures 3.1, 3.2, 3.3, 3.4, 3.5 and 3.6 are the double Gaussian fits for 
each of the three D° modes for either dataset, all in the Xp range 0.25 < Xp < 
0.35. 
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Figure 3.1: mD^-moo in GeV for B D*X with D^^K+n- in CLEO II 
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Figure 3.2: mD*-mDO in GeV for B D*X with D^^K+n- in CLEO II.5 
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Figure 3.3: mo^-mDO in GeV for B D*X with D^^K+n-n^ in CLEO II 
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Figure 3.4: mD*-mDO in GeV for B D*X with D^^K+Tr'Tr^ in CLEO II.5 
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Figure 3.5: mD*-mDO in GeV for B D*X with D^^K+n-n+n- in CLEO II 
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Figure 3.6: mD*-mDO in GeV for B D*X with D^^K+TT-n+n~ in CLEO II.5 



Shown in Table 3.3 are the uido and mD*-mDO double Gaussian data cuts for 
all three modes and for each dataset. 



Table 3.3: m^jo and mD*-'mr,o double Gaussian data cuts 







m£)0 




Mode 


Dataset 


in ± MeV 


in ± MeV 




CLEO II 


17.5 


1.15 




CLEO II.5 


15.0 


1.10 




CLEO II 


26.0 


1.50 




CLEO II.5 


27.0 


1.50 




CLEO II 


14.0 


1.30 




CLEO II.5 


12.0 


0.90 



3.6.3 Comparison with B D*X 

To verify the accuracy of our D* reconstruction, we compare our results with 
the PDG value Bi{B — D*X) = 22.7 % [5]. Since eMc varies with momentum, 
we need to compare our fully reconstructed sample with the expected number in 
the same momentum range. In D*~ p n and > D*~ p p n'^ signal Monte 
Carlo is ^ 50 % of the generated D*'s arc in the scaled momentum (Xp) range 
from 0.25 to 0.35. Approximately 32 % of Bi{B —>■ D*X) lies in this momentum 
range [40]. The expected D* yield is given by: 

(Number of charged/neutral B's) x 

(Br(B D* X) in X^ range 0.25 to 0.35) x 

(Br(D° mode)) x 

(Br(D* D^TTsoft) = 68.3 %) x 

{D* €mc) X 
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(0.32/0.50: efficiency correction). 

Tliis number is efficiency corrected for single and double Gaussian signal fits 
using > D*~ p n Monte Carlo. The corrected product is the expected yield 
that we compare to our respective fitted yield for each mode. 

The close agreement between the expected yield from Monte Carlo and the 
data result for each mode, as shown on Tables 3.4, 3.5, and 3.6, is a function of 
how well our Monte Carlo models D* decays from B mesons and of the accuracy 
of our reconstruction code. The double Gaussian fits yield better agreement than 
the single Gaussian fits. There is a significant drop in reconstruction efficiency for 
all three D° modes in the CLEO 11.5 data. This drop is due to the significantly 
reduced reconstruction efficiency of the soft pion in D* D^iTsoft in CLEO II. 5. 
Whereas in CLEO II the PTL drift chamber is used in the tracking algorithm, 
the SVX silicon detector which replaces it in the CLEO 11.5 data, in addition to 
having more material, is not used in the tracking algorithm. 



Table 3.4: D* yield: D^^K+tt- 



Dataset 


Single 
Gaussian 


Expected 
yield 


Found 


Double 
Gaussian 


Expected 
yield 


Found 


CLEO 11 
CLEO 11.5 


30.8 % 
16.3 % 


3,862.3 
3,933.1 


4,117.8 
4,455.2 


33.0 % 
20.7 % 


4,138,1 
4,994,9 


4,065.5 
5,083.6 
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Table 3.5: D* yield: D'^^K+tt-tt^ 





Single 






Double 








Dataset 


Gaussian 


Expected 


Found 


Gaussian 


Expected 


Found 


s/ 




<'A1C 


yield 




<'M(: 


yield 




VS + B 


CLEO II 


14.0 % 


6,371,3 


5,955.8 


15.3 % 


6,963.0 


6,856.6 


59.7 


CLEO II.5 


8.4 % 


7,356.0 


6,138.5 


8.9 % 


7,793.9 


7,531.3 


61.4 



Table 3.6: D* yield: D^^K+tt-tt+tt- 





Single 






Double 








Dataset 


Gaussian 


Expected 


Found 


Gaussian 


Expected 


Found 


s/ 




CMC 


yield 




CMC 


yield 




VS + B 


CLEO II 


14.5 % 


3,560.6 


3,637.6 


16.0 % 


3,928.9 


4,161.6 


48.3 


CLEO II.5 


6.5 % 


3,071.3 


3,100.5 


6.9 % 


3,260.3 


3,335.0 


47.7 



3.7 Antineutron Showers 

We need to define a set of criteria which allows us to select antineutrons with 
high accuracy without incurring a large loss in efficiency. The following character- 
istics, some limiting, some exploitable, apply to antineutron showers in the CLEO 
II electromagnetic calorimeter: 

1 . The antineutron shower yields an incomplete measurement of its energy and 
an accurate measurement of its direction. We use the antineutron shower 
energy to select candidates, but cannot use this energy when reconstructing 
the B candidate. We can use the well-measured shower direction. 
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2. Antineutrons frequently annihilate with matter in the calorimeter. Since 
antineutron annihilation showers have distinctive characteristics that enable 
us to separate them from other showers, we use these characteristics to select 
them. 

3. By baryon number conservation, if the event has an antineutron, it must have 
a corresponding baryon. In an exclusive reconstruction, as is ours, the selec- 
tion of an proton increases the probability of there being an antineutron in 
the event substantially due to baryon number conservation. In addition, once 
we have selected a D*, only shower candidates within a narrow momentum 
cone can be selected in the event. 

We refer the reader to other studies of baryon- antibaryon annihilation [43, 44]. 
In Table 3.7 we outline the types of showers encountered in a hadronic event. 

Table 3.7: Shower types and energy measured in calorimeter 



Particle 


Shower type 


Energy measured 


7r±, K^, p+ 

e^,7 

Kl 
neutron 
p and n 


minimum ionizing 
electromagnetic 
soft annihilation 
small electromagnetic 
medium to hard annihilation 


small fraction 
full measurement 

small fraction 
very small fraction 

small fraction 



We are limited by the absence of an independent antineutron sample that we 
can study to define our selection criteria. However, antiprotons also annihilate with 
nucleons in the calorimeter. Therefore, antiproton annihilation showers compose 
the shower sample in data and Monte Carlo which we use to define our antineutron 
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selection criteria as well as to gauge how well the Monte Carlo models antineutron 
annihilation showers. 

3.7.1 Shower Pcirameters 

We use the CCFC clustering package to select antineutron candidates. The 
CCFC shower package was not optimized to separate annihilation showers from 
other showers in the calorimeter. Nevertheless, we find that the parameters previ- 
ously optimized to include photons, such as E90E25, the list of nearby showers, 
and track-to-shower matching, are useful in excluding them as well. 

The shower parameters used are: 

1. E90E25 < cut 1, a pre-determined value: E90E25 has been defined in our 
discussion of 7r°'s. It is very close to 1 for photons, and, we find, farthest from 
1 for annihilation showers, overlapped 7r°'s- which are not merged, but have 
the two photon showers very near each other, and K^s. An energy dependent 
cut, cut 1, on E90E25, is applied to reject 99% of isolated photons. 

2. NNESH: List of nearby showers (which does not include the shower to which 
the list belongs). The area encompassed increases with energy. It is in this 
list that we find showers near the antineutron that are most likely to be 
hadronic split-offs. We call the shower associated with the list the main 
shower, and sum the energy of it and all the others in NNESH, we call a 
group. In D*- p n CLEO II signal MC, showers in the NNESH list 
are tagged to the parent shower in excess of 93 % of the time. 
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3. I {cos{9) I : The calorimeter can be divided by polar angle 9 into sections: good 
barrel, bad barrel, barrel/endcap overlap, bad endcap, good endcap, and near 
beampipe. We use the good barrel section throughout, which corresponds to 
\(cos(e)\ < 0.71. 

4. Match level: The track to shower match levels in CCFC's NTRSH array are: 

(a) Match level 1: shower center < 8 cm from the track projection. 

(b) Match level 2: not level 1, but > 1 member crystal 

(c) Match level 3: not level 1 nor level 2. 

We use match 3 showers for antineutron candidates and match 1 and 2 show- 
ers for antiprotons. 

3.7.2 Antiproton Showers in Data 

Lacking a sample of antineutrons to study in data which is independent of 
the sample we will use to measure D*" p n and D*~ p p tt"*" , we use 

antiproton annihilation showers which have been matched to a track. Antiprotons, 
like antineutrons, annihilate with matter in the calorimeter. The annihilation 
group of showers that results from this process has a distinctive signature: the 
antibaryon interaction is associated with a main shower, which has most of the 
energy of the annihilation, and the hadronic splitoffs from the interaction are 
associated with satellite showers near the point of interaction, each making a small 
contribution to the group energy. 

In Figure 3.7 we plot shower energy E^ain vs. signed momentum (PQCD) for 
protons and antiprotons, allowing all values of E90E25, in the CLEO II data from 
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a D* skim. Continuum and BB hadronic events are combined in this plot. The 
use of a D* skim is a result of convenience. 

The selection criteria used to generate Figure 3.7 is as follows: 

1. \PQCD{track)\ > 300 MeV. 

2. Lpr^Q^fyfi ^ 0.9. 

3. Track-to- shower match level 1. 

4. \crdE/dx,eiectron\ > 2.2 to rcjcct clcctron fakes. 

5. \(cos(9)\ < 0.71: good barrel 

6. All values of E90E25 allowed. 



In Figure 3.8 we use the same selection criteria as was used for Figure 3.7 with 
the exception of a rejection cut, E90E25 < cut 1, to suppress non-annihilation 
showers. The antiproton annihilation showers are typically those with E^ain > 
300 MeV and in the track momentum range 500 MeV to 900 MeV. Note that for 
these showers there is only a loose correlation between the reconstructed energy of 
the shower and the momentum of the antiproton that produced it. We can identify 
an annihilation shower, but have a very poor measurement of the momentum of 
the particle that produced it. 

The horizontal line at shower energy 200 MeV, sloping upwards at PQCD 1 
GeV, corresponds to minimum ionizing protons and antiprotons. The rest of this 
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ID IDB Symb Dote/Time Areo Mean R.M.S, 

782 1 1 990730/1023 2.8353E+04 - 1 .0707E-02 1-026 

0.5473 0.4285 

proton-onni. enrcc vs pqcd 



1.5 ■ ■, ^ ■ ,• 
1.0 
0.5 

0.0 r I r^-r 

-2.5 -1.25 1.25 2.5 

Figure 3.7: E^ain vs. PQCD for protons and antiprotons in CLEO II 

line in the lower momentum range corresponds to minimum ionizing pions. In the 
case of protons, before an E90E25 cut is applied, the diagonal line corresponds to 
protons captured in the calorimeter. 1.0 GeV momentum is the threshold for this 
capture. Neutrons, analogously to protons, do not annihilate, and, since we have 
neither momentum information nor distinguishable showers, the study of B decays 
with neutrons is beyond CLEO's capabihties. 

3.7.3 Antineutron Selection Criteria 

Using antiproton annihilation showers, we devise the antineutron selection cri- 
teria to be used when reconstructing > D*~ p n , as shown in Table 3.8. 
Although this selection criteria allows a shower to have no nearby daughters, in 
which case Emain = Egroup-i the majority of antineutron showers have at least one 
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ID IDB Symb Dote/Time Areo Mean R.M.S, 

782 1 1 990730/1023 l,3408E+04 -0,4700 0,9499 

0,8003 0.4182 



proton-onni. enrcc vs pqcd 



2.C 



1.^ 



1.C 



' ' 1 1 1,1 11 








■r. 




'it 
At 




•Mi 








i 





-2.5 -1.25 1.25 2.5 



Figure 3.8: Ej^ain vs. PQCD for protons and antiprotons in CLEO II 

nearby daughter in Monte Carlo, as well as in our exclusively reconstructed B' 
signal events. 

Table 3.8: Antineutron shower selection criteria. 



Track-to-shower match level 3 
E90E25 < cut 1 

\cOs{9shower)\ < 0.71 
Emain > 500 MeV 
Earoun > 800 MeV 



We test this selection criteria in a generic BB Monte Carlo before and after the 
requirement that there be a proton with \PQCD{track)\ > 300 MeV and Lp^otcm 
> 0.9 in the event. The generic Monte Carlo sample is discussed in Section 3.2. In 
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Table 3.9 we show the results before and after applying the antineutron selection 
criteria cuts for a sample before the proton requirement is applied. 

Table 3.9: Shower population in a generic BB Monte Carlo sample without proton 
requirement 



Shower 
type 


Annihilation cuts 
without Egroup cut 


Annihilation cuts 
with Egroup > 800 MeV 


7 from 7r° 


44.5 


46.2 


Kl 


33.0 


22.2 


P 


0.5 


0.6 




3.3 


1.9 


other (a;, etc) 


1.4 


1.3 


n 


17.3 


27.8 


TOTAL 


100.0 


100.0 



According to our generic Monte Carlo simulation, even before a proton require- 
ment is applied 27.8% of the annihilation-like showers in a BB event are antineu- 
trons. We next apply the proton requirement to a generic BB Monte Carlo sample. 
The results are shown in Table 3.10. 

Table 3.10: Shower population in a generic BB Monte Carlo sample with proton 
requirement 



Shower 


Annihilation cuts 


Annihilation cuts 


type 


without Egroup cut 


with Egroup > 800 MeV 


7 from 7r° 


16.3 


11.4 


Kl 


9.5 


3.9 


P 


2.7 


1.8 




2.0 


0.9 


other (a;, etc) 


0.7 


0.5 


n 


68.8 


81.5 


TOTAL 


100.0 


100.0 
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The above results are not surprising: use of baryon number conservation by 
selecting a proton increases the probability of an event having an antineutron from 
a few % to 50 %, while leaving all other annihilation like backgrounds at near the 
same level. Use of Egroup cut reduces Kl contribution by 1/3 in the Non-baryon 
sample and by 1/2 in the baryon sample. 

3.7.4 Antineutron Backgrounds 

In Figures 3.9 and 3.10 we compare the shower energy spectrum of each of the 
two major backgrounds to antineutrons-7 from tt^'s, and Kl-Io that of antineu- 
trons. All antineutron selection criteria cuts, including Egroup > 800 MeV, have 
been applied. Figure 3.9 is for case before proton requirement has been applied, 
and Figure 3.10 is for the case after proton requirement has been apphed. The 
sample is, again, generic BB Monte Carlo. The backgrounds that have not been 
included make up (4.3,3.0) % of the (before, after) proton requirement distribu- 
tions. In Figure 3.9 the solid distribution is for n's, the slanted lines distribution 
is for 7's from 7r°'s, and the white distribution is for K^s. 

In Figure 3.10 the white distribution is for n's, the solid distribution is for 7's 
from 7r°'s, and the slanted lines distribution is for K^s. The large contamination 
in the low range of the E^ain spectrum, as shown in Figure 3.9 and Figure 3.10 
does not affect D*^ p n since for this B meson mode the antineutron shower 

energy spectrum is dominantly in the range 1.0 < E^ain < 1-5 GeV. 
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Figure 3.10: Emain in GeV with proton requirement for n's, 7's from tt^'s, and 



CHAPTER 4 
MEASUREMENT OF 5° ^ D*-pPn+ 

Modes in S — > [D]NNX have not been previously measured. The reconstruc- 
tion efficiency multiphed by the number of B mesons in our dataset for D*~ 
p p 7r~^ and D*^ p n , assuming branching fractions similar to the ones mea- 
sured for S+ ^cPT^^ and B^ AcP7r+ tt" [29] , is at a level we can measure, 
which lends credibihty to a search for both modes. 

4.1 Monte Carlo Reliability 

The reconstruction process at CLEO, as well as at all other high energy de- 
tectors, involves a reconstruction efficiency eoata- For any given B decay, only a 
fraction of the number present in the data can be partially or fully reconstructed. 
If X number of B decays in our data are B Baryons decays, only eoata x ^ are 
measured. In the case of S — > pX, for example, there is a significant contamina- 
tion due to protons from beam gas interactions at the low end of the momentum 
spectrum which limits B pX to be carried out using antiprotons only, and then 
multiplying the result by 2. 

By tuning a wide range of parameters, the Monte Carlo simulation can be 
made to closely resemble the data, in which case the assumption eoata ~ ^mc is 
reasonable. €mc is rehable because several processes which are measured at CLEO 
allow considerable ease of simulation as well as large samples which can be reliable 
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separated from backgrounds. Processes such as e^e" — > e"^e~, e+e~ — > /I'^ii'', 
e^e~ 7r"'"7r~, e^e" —>■ 77, and e'^e~ — > 777 allow very accurate Monte Carlo 
modeling of individual particles, which in turn can be combined to model many B 
decays well. 

On this level of accuracy a second group of decays is studied to widen the scope 
of our simulation: Ks tt+tt", K~^K~, D* D^TTsoft with D^^K^ti~ , and 
A ^ pTT all have very low backgrounds. Many particle separation studies are based 
on these decays. An example is the use of a low background sample of A's to study 
antiproton annihilation showers, as will be discussed later in this chapter. 

4.2 Reconstruction Procedure 

B mesons are produced with 300 MeV momentum at CLEO, which implies 
that Eheam IS Very close to Eb- This momentum is small when compared with the 
B meson mass of 5.28 GeV. Since the beam energy at CLEO is measured with 
a 2 MeV resolution, we can constraint the mass of the B meson candidate to be 
equal to the beam energy Ebeam- We therefore use the beam contrained mass when 
reconstructing D*~ p p tt"*" : 

MbC = V^beam'-E-Pi" 

where Eheam is the beam energy, on average 5.29 GeV, and J2i Pi^ is the sum 
of daughter momenta squared. 

The energy difference between the beam energy and the energy of the recon- 
structed B candidate, defined by: 
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— -E'fceam " -^reconstructedB 

is centered at and has a Gaussian width similar to the Mbc Gaussian width, 
whereas background events are much more hkely to form a random distribution in 
Ae- Selecting B candidates using A^; cuts based on the Monte Carlo is useful in 
separating signal from background. As modelled by our Monte Carlo for 
D*~ p p n'^ is centered near zero and has a Gaussian width a < 100 MeV, which 
allows us to separate the signal from backgrounds that differ by a miss-measured 
extra pion. The methodology we use in measuring D*~ p p 7r+ is analogous 
to that used in previous exclusive reconstructions at CLEO [45, 46]. 

4.3 Monte Carlo Study 

In order to find the detector efficiency for the reconstruction of the modes we 
are searching for, we generate Monte Carlo in which one of the B mesons in the 
event is forced to decay to the mode we are reconstructing. We refer to this sample 
as signal Monte Carlo. In this signal Monte Carlo B mesons decay according to 
phase space. 

We also use a sample of generic Monte Carlo BB events on which we run 
the 5°— > D*~ p p TT"*" analysis code. The generic sample is composed of the B 
decays that have previously been measured as well as randomnly generated events 
assuming inclusive momentum distributions. There are no S — > [D]NNX modes 
in the generic Monte Carlo sample. We use the generic Monte Carlo sample to 
model our backgrounds due to modes of the D'^ we are not reconstructing. 
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Applying the particle selection criteria outlined in Chapter 3 for the > D*~ 
p p tt'^ decay daughters, we fit the AE distribution to a single Gaussian for signal 
and a 1st order Chebyshev polynomial for background. We show in Table 4.1 and 
Table 4.2 the result of these fits. 



Table 4.1: AE fit results for D*' p p tt+ signal MC in CLEO II 



Decay mode 


Central value (MeV) 


a (MeV) 


AE cut (MeV) 




0.6 


11.8 


± 35 




-2.3 


16.8 


± 50 




1.1 


9.8 


± 29 



Table 4.2: AE fit results for B^^ D* p p tt+ signal MC in CLEO II.5 



Decay mode 


Central value (MeV) 


a (MeV) 


cut (MeV) 




-1.0 


9.8 


± 29 




-3.1 


13.4 


± 40 




-1.6 


6.1 


± 18 



We define mode/dataset specific AE cuts which are applied to Mbc- The Mbc 
distribution is then fitted to a single Gaussian to determine each euc, as shown in 
Table 4.3 and Table 4.4. 

Table 4.3: Mbc distribution fit results for D*' p p n+ signal MC in CLEO 
II 



Decay mode 


eMC (%) 


a (MeV) 




8.69 ± 0.21 
4.48 ± 0.15 
3.89 ± 0.14 


2.63 
2.86 
2.65 
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Table 4.4: Mbc distribution fit results for D*' p p 7r+ signal MC in CLEO 
II.5 



Decay mode 


eMC (%) 


a (MeV) 




4.32 ± 0.22 
1.88 ± 0.14 
1.05 ± 0.08 


2.55 
2.60 
2.50 



The TTsoft in D* D^TTgoft for the CLEO II. 5 dataset is the most sensitive 
particle to errors in our Monte Carlo simulation. The replacement of the PTL 
detector by the SVX detector decreased by approximately a factor of 2 the recon- 
struction efficiency for the TTso/t's used here, and may have introduced a systematic 
uncertainty. 

4.4 Results in Data 

In Figure 4.1 and Figure 4.2 we show the AE vs Mbc distribution for > 
D*~ p p 7r+ in, respectively, CLEO II and CLEO II.5. In Figure 4.3 and Figure 
4.4 we show the Mbc distribution for > D*~ p p n'^ in, respectively, CLEO II 
and CLEO II. 5, after the mode dependent AE cuts have been applied. The mode- 
by-mode Mbc distributions are combined for each of Figure 4.3 and Figure 4.4, 
in which (dark) events are from continuum. The continuum background to > 
D*~ p p TT"*" is statistically insignificant. If this had not been the case, continuum 
subtraction would have been neccessary to subtract this background from the Mbc 
distribution. 
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Figure 4.1: AE vs Mbc distribution for D*' p p in CLEO II ON 

resonance data. 
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Figure 4.2: vs Mbc distribution for 5°^ D*' p p n+ in CLEO II.5 ON 
resonance data 
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Figure 4.3: Mbc (in GeV) for 5°^ D*- p p 7r+ in CLEO II 
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Figure 4.4: Mbc (m GeV) for 5°^ D*- p p 7r+ in CLEO II.5. 
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In Table 4.5 we attempt various fits to the combined CLEO II and CLEO II. 5 
Mbc distribution for D*~ p p tt~^ : 

1. Allowing Mbc and ctmbc float. 

2. Fixing Mbc = 5.28 GeV and allowing (Jmbc float. 

3. Allowing Mbc to float and fixing ctmsc according to expectations from signal 
MC. 

4. Fixing both Mbc and ctmsc prescribed in 2. and 3. above. 

The result of these fits are shown in Table 4.5. The variation in yield is small. 
Table 4.5: Results in the CLEO II/II.5 data of various fits to Mbc for D*- 



fitting options: 


fit results 


floating values 




Mbc (in MeV) 
a (in MeV) 
Fit yield 


5,280.7^^5 
30.52li|i 


fixed Mbc 




a (in MeV) 
Fit yield 


^•^U_o.53 

31.88t5'54 


fixed a 




Mbc (in MeV) 
a (in MeV) 
Fit yield 


5,280.5;jJ:^g 
2.65 
32.44l|68 


fixed Mbc and a 




Fit yield 
^(^0^ D*- ppn+) X 10-^ 


32.98^^:^^ 



B{B^^ D* p p 7r+ ) is calculated using: 
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B(B^^ D*- vvn-^ ) ^ 

O^n ly ppn ) ^2xB0B0xeMc)cLEOii+i2y<B'^B0xeMc)cLEOii.5 

The product {BB x ^mc)cleoii ~ product {BB x ^Mc)cLEOii.r>- 

The 

number of events for Mbc > 5.275 GeV is 21 for CLEO II, and 15 for CLEO II.5. 
Fitting each dataset separately would yield a measurable > D*~ p p n'^ ). 

In table 4.6 we quote the number of background events in the range 5.2 GeV 
< Mbc < 5.275 GeV, and the number of signal events for Mbc > 5.275 GeV. The 
D° mode with the largest background is >ir+7r~7r°. Signal region contains 
10 % background events. 



Table 4.6: B^^ D* p p : number of events found per mode 



Mode 


Background 


Signal 


Background 


Signal 




region 


region 


region 


region 




CLEO II 


CLEO II 


CLEO II.5 


CLEO II.5 




7 


6 





6 




14 


10 


15 


6 




4 


5 


3 


3 




5.225 5.250 5.275 

Beam constrained mass (GeV) 

Figure 4.5: Mbc for D*~ p p tt+ in data 



5.300 
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We find the A E distribution for S^— > D*~ p p tt^ , as shown in Figure 4.6 to 
be Gaussian, as expected from signal MC. We also find the momentum spectrum 
for TTs from the -D* decay, as shown in Figure 4.7, to agree with expectations from 
signal MC and to be very soft. In Figure 4.7 the solid distribution is CLEO II 
Monte Carlo and the dashed distribution is CLEO II/II.5 data. Only events with 
Mbc > 5.27 are plotted for both figures. 



MINUIT X Fit to Plot 
File: Generaled internally 



11-APR-2000 14:46 



X= 1.6for 14- Sd.o.f., 
Errors Parabolic 
Function 1: Chebysha/ Polynomial of Order 1 
NORM 280.17 ± 88.06 

CHEB01 0.36391 ± 0.3856 

Function 2: Gaussian (sIgma) 
AREA 21.387 ± 6.384 

MEAN -2.50699E-03 + 1.9945E-03 

SIGMA 6.57031 E-03 ± 1.5709E-03 

12 
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Figure 4.6: AE (in GeV) for B^^ D*~ p p Tr+ in CLEO II/II.5 data 
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500 




Figure 4.7: tt^ from D* momentum for > D* p p in data and Monte Carlo 
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4.5 Resonant Substructure 

In order to determine whether we are using the correct euc to measure B{B^^ 
D*~ p p TT'^ ), we need to know if we are measuring D*" p p 7r+ or a mode 
with intermediate particles decaying to the same combination of decay daughters. 
We therefore search for possible contributions to the resonant substructure of 
D*~ p p TT'^ . Any measurable resonant substructure would have to be subtracted 
from the inclusive measurement of B{B^^ D*~ p p tt"*" ), and the emc would 
likewise have to be adjusted. 

4.5.1 Two-body Decay and Possible Strong Resonances 
We search for two types of resonances: 

1. A heavy charmed baryon decaying strongly to p + D*~. No significant 
peaking is observed. 

2. A resonance of the virtual W decaying to ppn'^, > p p n'^ being a 
two body decay. No significant peaking is observed. 

We study the Mp^p^^+ spectrum in data to check if it is consistent with phase 
space or two-body decay. For the two-body sample we generate MC with a fictitious 
heavy particle that decays to ppTT"^ for >• D*~ p p 7r~^ . It has a mass of 2.6 
GeV and a width of 200 MeV. 

No conclusive evidence is found for a measurable contribution to B{B^^ D*~ 
p p 71^ ) from these possible contributions. The e^c for D*~ p p varies 

slightly depending on the resonance substructure. We allow a 5% systematic un- 
certainty to account for these variations. 



90 

4.5.2 A Baryon Contributions in the Form of S° D*-pA++ and S° 
D*-pA^ 

Only the non-resonant 5°— > D*~ p p ir'^ ^mc^ are shown in table 4.5.2. The 
CMC does not vary by more than 5 % for any of the three possible assumptions: 
D*- pp7r+ , ^ D*-pA++, and S° D*-pA^. Any contribution from 
the latter two would imply a lumping of events in the 1.0 to 1.3 GeV range for 
Mp^+ or Mp^+. The % makeup of Mp^+ and Mp^+ is nearly identical in Monte 
Carlo, and we quote only the former. There are 41 events in data in signal region. 

Tabic 4.7: % composition of generated (phase space MC) and events reconstructed 
(data) of Mp^+ and Mp^+ 



Mass 


Signal 


Data 




Data 




in GeV 


MC (%) 




% 




% 


1.0 - 1.1 


7.1 ± 2.1 





0.0 


2 


4.9 


1.1 - 1.2 


11.2 ± 0.8 


3 


7.3 


4 


9.8 


1.2 - 1.3 


11.9 ± 0.7 


8 


19.5 


7 


17.1 


1.3 - 1.4 


11.4 ± 0.6 


8 


19.5 


9 


22.0 


1.4- 1.5 


10.2 ± 0.6 


11 


26.8 


7 


17.1 


1.5 - 1.6 


9.2 ± 0.6 


3 


7.3 


6 


14.6 


1.6 - 1.7 


11.3 ± 0.6 


3 


7.3 


3 


7.3 


1.7- 1.8 


6.4 ± 0.6 


3 


7.3 


3 


7.3 



We place an upper limit of < 2 events for Mp^+, and < 3 events for Mp^+ 
as contributions from B^ — > D*~pA'^~^ , and B^ — > D*~pA^ respectively. These 
upper limits are educated guesses, and do not have any effect on the systematic 
uncertainty of our measurement of B{B^^ D*~ p p tt'^ ), since our quoted value 
includes the resonant substructure. Our quoted systematic uncertainty allows for 
the possibility of A baryon contributions to the resonant substructure of D*" 

P P TT'^ . 
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4.6 Backgrounds 

We use generic Monte Carlo, which is discussed in Section 3.2, to study possible 
backgrounds to D*~ p p n'^ . Several characteristics of > D*~ p p n'^ 

result in expected low backgrounds as well as an accurate Monte Carlo simulation: 

1. > D*~ p p TT"*" has a proton and an antiproton as decay daughters, the 
baryon-antibaryon constraint is strictly applied. 

2. I?*'s have low BB and continuum backgrounds. 

3. The high combined mass of the decay daughters is an added suppressant of 
continuum backgrounds. 

The AE cut suppresses modes of type > D*~ p p n'^ X, where X is any 
additional number of neutral and charged pions, but it does not suppress their 
combinatoric backgrounds. Since, as previously mentioned, B [D]NNX decays 
are not accounted for in the generic MC sample we study, it is reasonable to assume 
that if 75 % of the background in the range 5.2 GeV < Mbc < 5.275 GeV is not 
accounted for by the generic MC sample, the most likely source of this background 
is combinatoric background to B ^ [D]NNX decays. The small amount of non- 
B — > [D]NNX background is shown as a solid filled region in Figure 4.8, in which 
filled overlay on data is combinatoric non-B [D] NNX background from generic 
MC. 
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Figure 4.8: Mbc in data and generic Monte Carlo for D* p p tt"*" (in GeV) 
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4.7 Systematic Uncertainties 

In Table 4.8 we show the systematic uncertainties we consider to contribute 
significant errors to our measurement of > D*~ p p n'^ ). > D*~ p p n'^ 

has an average of 6.6 tracks. 

Table 4.8: Estimate of systematic uncertainties (in %) for 5°— > D*~ p p tt~^ 



Source 


Uncertainty (in %) 


branching fractions 


0.6 


D* branching fraction 


1.4 


D* reconstruction 


5.0 


Monte Carlo statistics 


5.0 


# of BB's 


2.0 


tracking (1%/track) 


6.6 


PRLEV proton ID 


8.0 


A signal contribution 


5.0 


A background contribution 


5.0 


Phase sj)ace axtsus two ])od\ 


5.0 


TOTAL 


15.5 % 



CHAPTER 5 
MEASUREMENT OF 5° ^ D*-PN 



As mentioned in Chapter 4, we have selected the B [D\NNX modes which 
are expected to have the lowest backgrounds as well as the highest reconstruction 
efficiency. When reconstructing > D*^ p n , however, we encounter a stumbling 
block not encountered in the reconstruction of > D*^ p p n'^ by virtue of the 
former mode having an antineutron-or a neutron in the case of its charge conjugate- 
as one of its decay daughters. Our inability to separate neutron showers from their 
backgrounds reduces the reconstruction efficiency for > D*~ p n by 50%. The 
antineutron reconstruction efficiency, however, is sufficiently large (in the range of 
30 to 40%). We do not have an accurate measurement of the antineutron energy, 
yet we do have a well measured direction of its shower, which we use to reconstruct 
> D*~ p n . The reconstruction method we use is analogous to that used in an 
unpublished reconstruction of i? ^ J/'^Kl by CLEO [47]. In Table 5.1 we outline 
our antineutron selection criteria which we derived in Chapter 3. 

Table 5.1: Antineutron shower selection criteria. 



Track-to-shower match level type 3 
E90E25 < cut 1 
|(cos(^)| < 0.71, or good barrel 

En^ain > 500 McV 

Egroup > 800 MeV 
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5.1 Reconstruction Procedure 



Instead of Mbc, we define m^o, in which we set AE = and assign the missing 
energy to the antineutron using the directional cosines of its shower. Since the 
measured shower energy for an antineutron, even after summing the energy in the 
hst of nearby showers, fails to match the total energy of the antineutron, we only 
use the electromagnetic shower energy as part of our selection criteria. 

The reconstruction steps for this mode are: 

1. En=Ebeam-ED''+p is assigned to the antineutron. This energy difference is 
the equivalent of setting AE = 0, or Ebeam = EreconstructedB- 



2. 3- momentum magnitude of antineutron candidate Pn—JE^ — m|, with nin 



= 0.9395 GeV. 

3. 3- momentum magnitude of antineutron times x,y,z directional cosines of 
shower energy are assigned to, respectively, x,y,z components of Pn. 

4. 4- momentum of B candidate p^o — Pn + Pd* + Pproton- 



As in the S'^— > D*~ p p n'^ reconstruction, we use Monte Carlo samples in 
which the B meson decays according to phase space. The resonant substructure 
of i)°^X+7r~7r° is same as in D*~ p p . 




5 



Mass of B candidate m^o = J E^^ 
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5.2 £)+ ^ p n in Monte Carlo 

In reconstructing D*' p n we can also be reconstructing B° — > Df D*" 
with Df p n. The B° Df D*' with Df p n contribution is part of our 
signal, yet it is not a mode within the resonant substructure of D*~ p n . 

It is a separate mode of the B meson with B{B^ D+ D*') = 9.6 ± 3.4 x lO'^ 
[5]. p n, a Feynman diagram of which is shown in Figure 5.1. Yet another 

contribution which can be in the Mp^n distribution is B^ — > D*^ D*~ with — > 
p n, which is a background. ^ D*+ D*~) = 2.0 ± 0.7 x 10"^ [5]. The effect 

of either on > D*" p n cannot be assesed accurately since we do not know 
B{D^ p n). If we assume B{D~^ — >• p n) is ~ 1%, the effect on >• D*~ p n 
can be significant enough to affect our measurement of > D*~ p n). 




Figure 5.1: A Feynman diagram for p n 
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The Mp+n distribution should not peak appreciably in > D*~ p n , and a 
narrow Gaussian peak with 5 MeV a in ^ Df D*~ with Df p n. The Df 
p n Monte Carlo peak for B'^ — > Df D*~ is shown in Figure 5.2. 




1.92 1.945 1.97 1.995 2.02 

Ds Mass (GeV) 

Figure 5.2: Mp+n in GeV. B° Df D*' with Df^pn 

In Figure 5.3 we show the smeared Mp^n distribution resulting from recon- 
structing B^^ D*~ p n in a signal Monte Carlo B^ D*^ D*~ with — > p 
n sample. Dfj — 1. The missing soft photon will cause this background 

sample to be considerably broader than the B° — > Df D*~ with p n signal 

contribution. 

The B^ D*+ D*~ contribution is a background to B^^ D*^ p n as well 
as B^ — > Df D*~ with Df p n. Our choice regarding how to deal with these 
contributions is to exclude both possible contributions (5° — > D*~ with 
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1.88 1.93 1.98 2.03 2.08 

Ds Mass (GeV) 



Figure 5.3: Mp^n (in GcV) from a reconstruction of D* p n in a signal 

Monte Carlo 5° ^ D*+ D*- with D+ ^ p n 

p n and 5° — > D*^ D*~ with — > p n) in quoting our final value for B{B^^ 
D*- pn). 

Both p n contributions can be excluded with 1.91 GeV < Mp+n < 2.04 

GeV. In Tables 5.2 and 5.3 we outline the ejuc's for 5°— > D*~ p n . These e^c's 
take into consideration that we are not reconstructing the conjugate mode-with a 
neutron. 
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200 




1.95 2.00 
Ds Mass (GeV) 



Figure 5.4: (white) is -> D+ D* 
with Df ^ pn. Mp+n (in GeV) 



with D-^ 



p n. (dashed) is 5° ^ D*+ D* 



300 



200 




Figure 5.5: (white) is B° Df D*' with Df 
with Df p n. uibo (in GeV) 



p n. (dashed) is S° ^ D* 
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Table 5.2: B^^ D*' p n signal MC smcs and widths in CLEO II. 



Mode 


% 


MeV 


both p n contributions excluded. 








7.47 ± 0.28 


3.10 




2.88 ± 0.17 


3.40 




3.33 ± 0.19 


3.11 


both p n contributions included, 






witliout any l)ackgrouii(l correctioii. 








8.12 ± 0.32 


3.14 




3.09 ± 0.19 


3.00 




3.52 ± 0.20 


3.05 



Table 5.3: D*' p n signal MC cmc's and widths in CLEO II.5 



Mode 



% 



MeV 



Both D+ 



p n contributions excluded. 



K+n 

n 







4.51 ± 0.17 

1.60 ± 0.13 
1.68 ± 0.14 



2.95 
3.23 
2.79 



Both p n contributions included, 

without any background correction. 



4.99 ± 0.18 

1.82 ± 0.14 

1.83 ± 0.15 



2.92 
3.35 
2.78 
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5.3 Results in Data 

We look for > D*' p n with or without the mass region 1.91 GeV < Mp+n 
< 2.04 GeV. We compare in Figure 5.6 these two cases in CLEO II/II.5. 



16 



I — ^ 
5.200 



5.225 



5.250 
B mass (GeV) 



5.275 



Figure 5.6: (white) is inclusive of all contributions, (dashed) is after exclusion of 
both p n contributions 



The 5"— > D*~ p n ttibo distribution is less statistically significant without the 
Df p n contributions. However, the cost of lost events for the purposes of our 
measurement is warranted for the following reasons: 

1. The signal events we are left with have an insignificant probability of being 
something other than D*~ p n . 

2. We properly account for the lost events in our eMc's- 

In Figure 5.7 we show mso for ON resonance and OFF resonance for the case 
we use to quote B(-B°— > D*' p n): without both contributions. 
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Figure 5.7: (white) is ON resonance, (solid) is OFF resonance 

The number of OFF resonance events is not significant enough to affect our 
measurement. We choose not to subtract these from the ON resonance distribution. 
In Figure 5.8 and Figure 5.9 we show itibo distributions by dataset. 



File: *bnnd>!_apr4_20.his 

ID IDB Symb Date/Time 
539 5 1 000411/1602 



5.200 5.225 5.250 5.275 5.300 

B mass (GeV) in CLEO II 
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Figure 5.8: Mgo for D*' pn m CLEO II ON resonance data 



File; *bnndx_apr3_25.his 

ID IDB Symb DateAlme 
539 10 1 000411/1604 



E 4 



5.200 



n nnl 



LA 



5.225 5.250 5.275 

B mass (GeV) in CLEO ii.5 



5.300 



Figure 5.9: Mgo for D* p n in CLEO II. 5 ON resonance data 



104 

5.4 — > p n in Data 

In Figure 5.10 we plot Mp+^ for the region where we expect to see Df p n 
as a small width signal for B° Df D*~ , and as a shallow background for 

D*~ . The horizontal lines demarcate the mass region 1.91 GeV < Mpj^n 
< 2.04 GeV, which we are excluding. Events with meo > 5.275 GeV region are 
shown in Figure 5.10. The vertical hues demarcate the mass region 1.91 GeV < 
Mp+ji < 2.04 GeV, which we are excluding. Not shown are events with Mp+^ > 
2.27 GeV. We expect to find ^ 10% of > D*~ p n signal events in this region, 
which translates to 3 events. 



E 




1.87 1.97 2.07 2.17 2.27 

mass of proton + antineutron in CLEO II/II.5 



Figure 5.10: Mp+n for 5"-^ D*- p n (in GeV) in CLEO II/II.5 

We find 8 events where a Df —>■ p n signal is expected. This sample is not 
sufficiently significant to label these events as signal. The scatter plot in Figure 
5.11 contains more information than Figure 5.10: there are other regions in this 
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plot that are equally statistically significant to the region where we expect to find 
p n. If we are to claim a signal in one, we should be able to claim one in 
any of these others. 



2.67 



S 2.47 



+ 
c 
o 

Q 
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1.87 I — ' — ' — ' — ' — ' — ' — ' — ' — ' — ' — ' — ' — ' — ' — ' — ' — ' — ' — ' — I 
5.200 5.225 5.250 5.275 5.300 

B mass (in GeV) in CLEO II/II.5 

Figure 5.11: Mp+^ vs. m^o (both in GeV) 

We can, however, estimate an upper limit for Df pn. The most conservative 
estimate is to call all 8 events in the expected region signal, which leads to a value 
of 11.5 for a 90% confidence level measurement. 11.5 signal events correspond to 
B{Df ^pn) < 5.4 % using the PDG value for B{B^ Df D*-). 
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5.5 B(B^^ D*- p n ) Measurements. 

We use the same equation to calculate B{B^^ D*~ p n ) as we did for B^- 
D*~ p p 7r+ : 



B{B^^ D*- p p 7r+ ) = 



FittedYield 



^Mc)cLEOII.l 



As we mentioned when calculating B{B'^^ D* ppn'^), the product {BB x eMc) 
is almost the same for CLEO II and CLEO II. 5. The number of events for this 
mode for Mbc > 5.275 GeV is 10 for CLEO II, and 16 for CLEO II.5. 

The combined uibo distribution can be fitted to yield results that are consistent 
with Monte Carlo expectations. 




5.200 



5.225 5.250 5.275 

B mass (GeV) 



5.300 



Figure 5.12: Inclusive m^o in data for > D* p n (in GeV) 
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5.200 5.225 5.250 5.275 5.300 

B mass (GeV) 



Figure 5.13: uibo in data for > D* p n in GeV excluding both p n 

contributions 



Table 5.4: Measurements of > D* p n branching fractions 



Mode 


Yield 


Branching Fraction 


> D*" p n inclusive 
> D*~ p n without both contributions 


32.941«:^^ 


14.911^ X 10-^ 

ii.9li:i X 10-' 



The difference in branching fractions does not imply that what has been ex- 
cluded is only p n. We do not attempt to derive an upper limit by this 
method. As in the case of our measurement of D*~ p p 7r~^ , we find the 
D*~ p n signal events are consistent with phase space decay. We look for a 
resonance of type {cud) + {dd)^^^^^ — > n + D*~ , and fail to find it. 

The fits are quoted in Table 5.5. mso is fitted to a single Gaussian for signal, 
and an Argus function for background. The offset is fixed at 0.0025, and the £^6eam 
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is fixed at 5.29 GeV in all the fits. As in the case of > D*~ p p n'^ , in table 
5.5, leaving some parameters fixed, but not others allows us to test the validity of 
our results. The excluded region is 1.91 GeV < Mp^^ < 2.04 GeV. 

Table 5.5: Results of various fits to tubo for B^—>^ D*~ p n with or without the 

p n contributions 



Fitting options 


No p n 


Inclusive 


floating values 






niBO (MeV) 
a (MeV) 
Fit Yield 


5,280.3l[j:^^ 

2.59in? 
23.57ltJ^ 


5,280.51[J:^^ 
2.43l°i? 
31.48l|i 


fixed mfio = 5.28 GeV 






a (MeV) 
Fit Yield 


9 fi9+U.4y 
^•"^-0.42 

23.67ll:i 


9 c;c;+0.47 
^.OO_0.40 

31.97l|i 


fixed a = 3.07 MeV 






m^o (in MeV) 
Fit Yield 


5,280.3ll]:r7 
23.97111^ 


5,280.31^0 
32.75lij? 


fixed m^o and a 






Fit Yield 


24.021^:^^ 


32.941^:^^ 


B{B^^ D*- pn) X 10"^: 






From raw yield 


ii.9ii| 


14.911? 


With MC Correction Factor 


14.8l|f 


I8.3II3' 



In the following section we discuss the correction factor we apply to the raw 
B{B^^ D*~ p n ) due to a discrepancy between data and Monte Carlo in the ep for 
antiproton annihilation showers. The Correction Factor we use increases the raw 
B{B'^^ D*~ p n ) quoted in Table 5.5 by a factor of 1.21 with an added systematic 
uncertainty of 4%. 
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5.6 Correction Factor 

The Monte Carlo we use has not been optimized to model nucleon-antinucleon 
annihilation. We find a discrepancy between the reconstruction efficiency for data 
and Monte Carlo for antiprotons. We consider the discrepancy credible and change 
the antineutron efficiency using a Correction Factor. We assume that the Monte 
Carlo fails to model antineutrons by the same amount as it does for antiprotons, 
even though the quark content of antiprotons and antineutrons is different. 

We define eannihUation as the efficiency for an antibaryon shower to pass the 
annihilation shower cuts outlined in Table 3.8. This quantity will be different for 
antiprotons and antineutrons. In the former case we choose a sample which is 
inevitably biased by the momentum range in which we can separate antiprotons 
from other charged tracks. In the latter case we use a sample which has shower 
backgrounds. 

We define a correction factor to account for the discrepancy between data and 
Monte Carlo for annihilation showers as: 



Correction Factor = 

^annihilation, Data 

This Correction Factor is found for antiproton momentum bins in a range 
consistent with the momentum range of antineutrons in > D*~ p n as found 
from the energy assigned to the antineutron candidate shower after assuming the 
antineutron mass. The number of antineutrons in each of these momentum bins 
is used to weigh the contribution of each bin to the Correction Factor, which is 
defined for the entire range. 
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5.7 Use of a A Sample 

The selection of a clean sample of antiprotons is achieved by reconstructing A's. 
In the decay A Ptt"*", due to the large difference in mass between the p and the 
TT"*", the kinematics of the decay allow for a clear separation of the p tracks from 
the TT"^ tracks. Furthermore, the single Gaussian width of the distribution is 1 
MeV and is affected by backgrounds insignificantly 

Antiprotons are selected with the same criteria used in the reconstruction of 
> D*~ p p 7r+ except for the use of Lproton > 0.1. A's are from the KNVF 
package [42]. The A selection criteria is as shown in Table 5.6. 



Table 5.6: A selection criteria 



Both tracks pass TRKMNG flag tng(track) > 

of two tracks to form a vertex < 30 
for A to point back to the interaction region < 30 
Significance of the two dimensional flight distance < 3(7 
Flight distance > 0.005 meter 



We count the number of A's in the signal region after applying a double Gaus- 
sian fit with fixed parameters as shown in Table 5.7. 

Table 5.7: A double Gaussian fixed parameters 



Mass 


1.1158 GeV 


a' 


0.00162 


Arca2 / areal 


0.46 


'^MEAN 





ay a' 


0.425 



Ill 

In Table 5.8 we find the % of antiproton showers in a A sample which pass 
the annihilation shower cuts. In Table 5.8 we find eannihUation for MC and Data A 
samples. The samples are a combination of ON and OFF resonance. The ON and 
OFF resonance eannihUationS were found to yield the same results. 

Table 5.8: e annihilation'- {%) of match 1 or 2 antiprotons in A's passing annihilation 
shower selection cuts for six momentum spectra 





CLEO II 


CLEO II 


CLEO II.5 


CLEO II.5 


Momentum range 


MC 


Data 


MC 


Data 




% 


% 


% 


% 


300-500 MeV 


28.3 ± 1.3 


13.2 ± 0.3 


26.5 ± 1.0 


13.3 ± 0.2 


500-700 MeV 


65.1 ± 1.7 


52.7 ± 0.9 


65.0 ± 2.6 


52.2 ± 0.6 


700-900 MeV 


70.8 ± 1.8 


60.6 ± 1.0 


71.5 ± 2.8 


60.9 ± 0.7 


900-1100 MeV 


65.3 ± 1.7 


52.9 ± 0.9 


63.4 ± 2.5 


50.7 ± 0.6 


1100-1300 MeV 


50.6 ± 2.0 


41.9 ± 1.0 


52.2 ± 2.1 


41.4 ± 0.7 


1300-1700 MeV 


47.4 ± 1.9 


44.0 ± 1.0 


46.3 ± 1.8 


41.0 ± 0.7 


1700-2100 MeV 


56.0 =b 2.0 


54.0 ± 1.0 


55.1 ± 1.9 


53.2 ± 0.7 



The Correction Factor is weighed according to the fraction of antincutrons in 
a momentum range as produced in Monte Carlo. The solid distribution in Figure 
5.14 is data and the dashed distribution is signal Monte Carlo. 
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File: Generated internally 

ID IDB Symb Date/Time Area Mean R.M.S. 

539 11 1 000417/1404 26.00 1.108 0.2619 

539 4 2 000417/1407 26.01 1.077 0.3316 




Figure 5.14: Antineutron momentum distribution (in GeV) 
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In Table 5.9 we calculate the Correction Factor as previously defined, as well 
as the fraction of antineutrons found in signal Monte Carlo for each momentum 
range. 

Table 5.9: Correction Factor (C.F.) as a function oipproton and % of antineutrons 
found in signal Monte Carlo for each momentum range in CLEO II/II.5 



Momentum range 


% of n in MC 


C.F. 


300-500 MeV 


9.1 


2.07 ± 0.11 


500-700 MeV 


10.0 


1.25 ± 0.04 


700-900 MeV 


14.7 


1.17 ± 0.04 


900-1100 MeV 


15.2 


1.24 ± 0.04 


1100-1300 MeV 


16.3 


1.23 ± 0.06 


1300-1700 MeV 


20.7 


1.10 ± 0.05 


1700-2100 MeV 


15.7 


1.04 ± 0.05 


300-1700 MeV 


100.0 


1.21 ± 0.05 



We will use 1.21 as the Correction Factor by which to increase B{B^^ D*~ p 
n ), which is an equally weighed average of CLEO II and CLEO II. 5, and add to 
our list of systematic errors a 4% contribution due to this correction. 
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5.7.1 Backgrounds in Signal and Generic Monte Ceirlo 

To check that our antineutron selection criteria does not contaminate our sam- 
ple with showers from decays similar to D*~ p n ,we run our > D*~ p n , 
with D^^K^t:~ reconstruction code on signal Monte Carlo for D** p n with 
equal numbers of £11(2420)° and L'2*(2460)°, B^ D*-pn7r^, and B^ D*-pA^ 
with A° — > mr^. The result for m^o are overlayed in solid on our data results, in 
white, for > D*~ p n m Figure 5.15. 



300 - 



200 - 



100 - 



JO 



Figure 5.15: uibo in GeV for > D* p n selection criteria applied to selected 
signal Monte Carlo background modes 



We also run the > D*~ p n code for all three modes on our generic 
Monte Carlo sample, as shown in Figure 5.16. The generic Monte Carlo sample is 
discussed in Section 3.2. Scaling the generic Monte Carlo sample to the size of our 
dataset, we find it to account for 40% of the background. 
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Within the generic MC sample f^i 25% is combinatoric background events from 
each of D^^K'^n'Tr'^n' and D'^^K'^n' . The remaining events, 50%, are from 
2)°— ^i^+Tr"?!^. Since B — > D*^NNX decays are not included in the generic MC, 
these events are from random combinations of a S — > Baryons decay and a S — > 
D* X decay. 




5.200 5.225 5.250 5.275 5.300 



Figure 5.16: m^o for > D* pn m GeV. (white) is data, (solid) is generic MC 

We are unable to account for 60% of the background found in data in the m^o 
distribution for > D*~ p n . However, this background is small, and the m^o 
distribution we find is consistent with a signal at the nominal B mass. 

5.8 Antineutron Directional Cosine Resolution 

We compare the generated and reconstructed directional cosines for the an- 
tineutron candidates in signal Monte Carlo. We find the median error to be 37.8 
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milliradians. Since a single 5-cm calorimeter cell corresponds to ~ 50 milliradians, 
the systematic error in our measurement of the directional cosines of the antineu- 
tron candidate shower is insignificant. 

5.9 D*+pn. 

Neutrons do not annihilate in the calorimeter. Since our selection criteria for 
antineutron showers rejects the vast majority of neutrons, we derive a new set of 
cuts specific to neutrons. The selection criteria we used is shown in Figure 5.10. 
The most important difference with antineutrons is the selection of low energy 
photon-like showers. 



Table 5.10: Neutron shower selection criteria 



Track-to-shower match level type 3 
E90E25 > cut 1 
\{cos{6)\ < 0.71, or good barrel 

Emain < 500 McV 



We are unable to define selection criteria that allow us to separate neutrons 
from soft photons. This background overwhelms the signal, making unatainable 
the reconstruction of 5° D*~^pn. We plot E^ain for neutrons (solid distribution) 
and antineutrons (dashed distribution) in Figure 5.17. 
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File: ■/diskl/oir/bnndVneutfon.slgtioltnc.his 

ID IDB Symb Dote/Time Area Mean R.M.S. 

300 1 1 991021/1048 712.0 0.1664 0.1451 

300 2 2 991021/1048 653.0 '.CB8 0.4571 




Figure 5.17: Emain (in GeV) for neutrons and antineutrons in > D* p n signal 
Monte Carlo 
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5.10 Systematic Uncertainties 

In Table 5.11 we show the systematic uncertainties we consider to contribute 
significant errors to our measurement of > D*' pn). > D*~ p n has an 

average of 4.6 tracks. 

Table 5.11: Estimate of systematic uncertainties (in %) for > D*~ p n 



Source 


Uncertainty (in %) 


branching fractions 


0.6 


D* branching fraction 


1.4 


D* reconstruction 


5.0 


Monte Carlo statistics 


5.0 


# of BB's 


2.0 


tracking (1%/track) 


4.6 


PRLEV proton ID 


4.0 


A background contribution 


5.0 


Phase space versus two body 


3.0 


antineutron miss-ID 


15.0 


TOTAL 


18.8 % 



CHAPTER 6 
CONCLUSION 



This work comprises the second successful exclusive reconstruction oi B ^ 
Baryons modes, after the exclusive reconstructions of S+ — > AcPTr"*" and S° — > 
AcPTT"*" TT~ carried out by CLEO in 1997 [29]. In Table 6.1 we summarize all the 
exclusive B — > Baryons measurements to date. We rank the modes from largest 
to smallest branching fraction central value. The first error is statistical and the 
second is systematic. The Ac modes have a third uncertainty, which is systematic 
and is due to the error in the determination of B(A+ pK~ 7r~^). 

Table 6.1: Exclusive measurements in B ^ Baryons to date 



Mode 


Branching Fraction (xlO ^) 


5°^ D*- p n 

^ AcP7r+ TT" 

B+ AcP7r+ 


14.51^:^ ± 2.7 
13.3l4:2±3.1 ±2.1 

6.61^:1 ± 1.0 
6.2l|^ ± 1.1 ± 1.0 



As shown in Table 6.1, the size of the branching fractions we are measuring 
in this work are of the same order of magnitude as those previously measured in 
decay modes including a Ac. We are unable to extrapolate the exclusive results 
to compare the relative magnitude of the inclusive modes with or without a Ac. 
For instance, the available phase space for the production of more light mesons 
is larger in the case of modes with a Ac than in the case without, which may 
cause the B AcX inclusive rate to be larger than the B — > [D]NNX inclusive 
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rate. However, we think it reasonable to assume that the measurements we have 
at our disposal to date are representative of the inclusive modes. We infer that 
B —>■ [D]NNX modes contribute significantly to the total B — > Baryons decay 
rate. 

6.1 5 — > Baryons phenomenology 

The B Baryons theoretical models attempted to date [48, 49] rely on the 
assumption that the B meson decays to two decay daughters, which is the two- 
body assumption. Unlike two-body B meson modes, most B Baryons modes, 
including all four that have been measured to date, have at least three decay 
daughters. Most B Baryons modes also have a varied resonant substructure. 
The assumption that the hadronization process takes place at a late stage in the 
decay, which allows for a substantial simplification of the equations due to the 
suppressed dependence on the exchange of gluons and light quarks, is justified in 
hard two-body decays. However, in many-body decays there can be re-scattering 
of the hadronizing quarks and multiple exchanges of soft gluons and virtual quarks. 
Notwithstanding the calculational difficulties of many-body B Baryons decays 
such as the ones we measure here, progress continues to be made using HQET to 
explain B meson decay. The focus has been on decays for which HQET can be used 
to extract the most information, and which can be used to search for CP violation 
[14, 15]. Future attempts to theoretically explain B — > Baryons are needed. 
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6.2 Possible Future B Baryons Modes at CLEO 

In CLEO II reconstruction from B decays is plagued by high backgrounds 
due to poor separation between kaons and pions. The substantial improvements in 
charged particle separation (pion-kaon-proton) in the CLEO III data should allow 
for the successful reconstruction of S+ — > D^pn and B^ — > D^ppn'^. 

The ability to reconstruct antineutron modes should allow for the reconstruc- 
tion of modes of type B — > (Ac, Sc)nX. Modes of type B — > (Ac*, Ec*)pX, and 
B — > ScAc, which have low reconstruction efficiencies in CLEO II/II.5, are also 
worth pursuing. 

B — > Baryons for b ^ sg, and b ^ u modes will also be of interest and likely 
within range. If previous measurements of the meson modes are used as a guide, 
the modes S+ — > Ap, B^ — > pnTr~, and S° — > ppTr'^Tr~ are worth pursuing. 

Semileptonic b ^ c B ^ Baryons decays, such as B'^ — > Ac~pl~^i'i, with I = e, /i, 
due to the low reconstruction efficiency at CLEO for low momentum leptons, have 
not been measured. The substantially larger dataset expected in CLEO III may 
compensate for the low reconstruction efficiency. A mode such as B^ — >■ ppl^vi 
should be feasible. This latter mode can be used as an auxiliary in the measurement 

of iKbl- 

6.3 Significance of Results 

We have found the first evidence of decay modes of the B^ of the type B — > 
DNN-K. We measure the branching fractions B (S°^ D*~ p p 7r+) = (6.5li;2 ± 
1.0) X 10-^ and i3( D*- pn ) = (14.811^ ± 3.0) x 10"^ These measurements 
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indicate the fraction of baryonic decays of B mesons that do not proceed via A+ 
may be of approximately the same magnitude as those that do. 
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